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Abstract: We consider large U{N)^ thermal M = 2 quiver gauge theories on x . 
We obtain a phase diagram of the theory with i?-symmetry chemical potentials, separating 
a low-temperature/high-chemical potential region from a high-temperature/low-chemical 
potential region. In close analogy with the = 4 SYM case, the free energy is of order 
C(l) in the low-temperature region and of order 0{N'^M) in the high-temperature phase. 
We conclude that the M = 2 theory undergoes a first order Hagedorn phase transition at 
the curve in the phase diagram separating these two regions. We observe that in the region 
of zero temperature and critical chemical potential the Hilbert space of gauge invariant 
operators truncates to smaller subsectors. We compute a 1-loop effective potential with 
non-zero VEV's for the scalar fields in a sector where the VEV's are homogeneous and 
mutually commuting. At low temperatures the eigenvalues of these VEV's are distributed 
uniformly over an jl^u which we interpret as the emergence of the /I^m factor of 
the holographically dual geometry AdS^ x /I^m. Above the Hagedorn transition the 
eigenvalue distribution of the Polyakov loop opens a gap, resulting in the collapse of the 
joint eigenvalue distribution from /Zim x into I'Lm- 
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1. Introduction 



The phase structure of large N U{N) gauge theories at finite temperature is in itself a very 
rich and interesting subject that may provide qualitative insight into the phase structure of 
QCD. Even more so, the AdS/CFT correspondence ||l|, ^ has provided a general frame- 
work for translating results obtained in weakly coupled thermal gauge theory into results 
about the finite temperature behavior of the physics of black holes and stringy geometry 
at strong coupling. One such connection was suggested by Witten ^ who argued that the 
Hawking-Page phase transition |p between thermal AdS^ and the large AdS^ Schwarzschild 
black hole should have a holographic dual description as a confinement / deconfinement tran- 
sition in the dual thermal field theory defined on the conformal boundary x of thermal 
AdS5. 

A general framework for studying large N U{N) gauge theories on at finite tem- 
perature was given in ||6|. In particular, this considered = 4 U{N) SYM theory which 
was also independently studied in Furthermore, for the AA = 4 case the analysis was 
extended in ||8|, |9|] to include chemical potentials conjugate to the i2-charges. In this way a 
phase diagram of the theory as a function of both temperature and chemical potentials was 
obtained. As one application of the phase diagram, in the observation was made that in 
regions of small temperature and critical chemical potential A^ = 4 SYM theory reduces to 
quantum mechanical subsectors, including the XXX1/2 Heisenberg spin chain. ^ 

Again for AA = 4 SYM theory, the framework of |^ was generalized in a different 
direction in by allowing non-zero VEV's for the scalar fields. There a one-loop effective 
potential for the theory at finite temperature on 5^ at weak 't Hooft coupling was computed 
under the assumption that the VEV's of the scalar fields are constant and diagonal ma- 
trices.^ The potential was used there to study the manifestation of the Gregory-Laflamme 
instability^ of the small AdS^ black hole from the weakly coupled gauge theory point of 



view. The solutions to the equations of motion obtained from the effective potential of 



were given in |13| in terms of a joint eigenvalue distribution of the Polyakov loop and the 
scalar VEV's. Within the sector of constant and commuting scalar VEV's it was found that 
the topology of the eigenvalue distribution of these VEV's undergoes a phase transition 
X ^ at the Hagedorn temperature. The authors interpreted the eigenvalue dis- 
tribution of the scalar VEV's as the emergence of the factor of the holographically dual 
thermal AdS^ x geometry. It should be noted that, while the truncation to commuting 
matrices is consistent, this sector will not describe the absolute minima of the action |14|. 



For this reason the observed phase transitions in the commuting saddles studied in Ref. ||T3| 
are not transitions in the full gauge theory. 

The discovery that the eigenvalues of scalar VEV's reconstruct the dual spacetime ge- 
ometry was originally made by Berenstein et al. in |15, [l^, |l8| by setting up matrix 



models for the various sectors of BPS operators in the chiral ring. In particular the model 



^Recently other decoupling limits have been found in near-critical regions by extending this analysis to 
include the chemical potentials conjugate to the angular momenta on S'"^ p^. 

^This potential was computed earlier in for the special case of zero Polyakov loop eigenvalues. 
^See ||l2|| for a recent review of the Gregory-Laflamme instability. 
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for 1/8 BPS operators was developed in |17| where the dynamics was shown to reduce to 



that of the eigenvalues of three commuting Hermitian matrices Z, X, Y plus two fermionic 
matrices^ Wa- The quantum mechanical Hamiltonian for the eigenvalues involves an at- 
tractive harmonic oscillator part and a repulsive Vandermonde type part. These forces are 
balanced when the eigenvalues are localized to a hypersurface in which is taken to be an 
due to the 50(6) invariance of the quantum Hamiltonian. This was identified with 
the factor of the holographically dual geometry AdS^ x S^. 

The purpose of this paper is to investigate the phase structure of = 2 U{N)'^'^ quiver 
gauge theories at finite temperature.^ Defined at zero temperature and on a flat spacetime, 



these gauge theories are J\f = 2 supersymmetric and conformally invariant 25 1. We carry 
out the investigation of the phase structure in two directions. First, we consider the case 
of non-zero i2-symmetry chemical potentials. One interesting question here is whether the 
high-temperature phase admits several solutions. A further point of interest is to examine 
whether one can uncover information about closed subsectors of the as yet not completely 
settled underlying spin chain of = 2 quiver gauge theory by studying the near-critical 
chemical potential and low temperature regions of the (T, n) phase diagram of the theory 
as done for 7^ = 4 SYM theory |, |2|. 

Another question of interest is to what extent the 5^ eigenvalue distribution of the 



AA = 4 SYM scalar VEV's found in |17| and |13| can be interpreted as the emergence of 
the factor of the dual string theory geometry AdS^ x 5^. To examine this question, we 
make use of the fact that M = 2 quiver gauge theory can be realized as a Zm projection 
of = 4 SYM theory. The holographically dual spacetime of the M = 2 theory is thus 
AdS^ X jTLu where Za/ only acts on the factor. If the above interpretation of emergent 
spacetime is correct, we should then expect to find an jTLu eigenvalue distribution for 
the VEV's of the scalar fields of AA = 2 quiver gauge theory. This has been studied via 
counting of BPS operators in |27, 29|. Our approach to the problem is complementary 
in that it is valid for weak 't Hooft coupling, and it is valid for all temperatures in the range 
< Ti? < unlike 1^, |8|, Hi which is only vahd for T = 0. In parallel with Ref. 



we restrict to the sector of constant and commuting scalar VEV's. Whereas this enables us 
to study phase transitions in the eigenvalue distributions, revealing interesting dynamics, it 
does not necessarily reflect the full phase structure. However, we find it enlightening to see 
how the geometry of the dual AdS spacetime is mirrored in the structure of the quantum 
effective action computed in this sector. 

The outline and summary of the results in this paper are as follows. In Section 2 we 
give an introduction to AA = 2 quiver gauge theory on x with chemical potentials 
conjugate to the i?-charges. In Section 3 we evaluate the quantum effective action oi M = 2 
quiver gauge theory with non-zero i?-symmetry chemical potentials and zero scalar VEV's 
in the gym limit and express it in terms of single-particle partition functions. We use 
the effective action to construct a matrix model for AA = 2 quiver gauge theory on X S^. 
The model turns out to be an M-matrix model with adjoint and bifundamental potentials. 



However, throughout the analysis of the dynamics in the fermionic matrices Wa are disregarded. 
See also Refs. ||l^, |2^, ^ for related work on other supersymmetric gauge theories. 
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In Section 4 we study the saddle points of the matrix model as functions of tempera- 
ture and chemical potential and thereby examine the phase structure of the model. In the 
low-temperature phase we find a saddle point corresponding to a uniform distribution of 
the eigenvalues of the Polyakov loop^. In this phase the free energy is 0{1) with respect 
to N. This behavior of the free energy suggests that the model in this phase describes a 
non-interacting gas of color singlet states, and the phase is therefore labelled "confining". 
This saddle point is observed to become unstable when the temperature is raised above a 
certain threshold temperature (which depends on the chemical potential). The model then 
enters a new phase in which the free energy scales as N'^M as N ^ oo. This phase is 
thus interpreted as describing a non-interacting plasma of color non-singlet states and is 
labelled "deconfined". The "deconfinement" transition is of first order and is identified with a 
Hagedorn phase transition. The condition of stability of the low-temperature saddle point is 
translated into a phase diagram of the gauge theory as a function of both temperature and 
chemical potentials. We subsequently study the phase diagram in regions of small temper- 
ature and critical chemical potential. We observe that the Hilbert space of gauge invariant 
operators truncates to the SU{2) subsector when the chemical potential corresponding to 
the SU (2)/j factor of the i?-symmetry group SU{2)r x U{1)r is turned on, whereas when 
both chemical potentials are turned on and set equal, it truncates to a larger subsector that 
corresponds to an orbifolded version of the SU {2\3) sector found in = 4 SYM theory. 

In Section 5 we develop a matrix model for M = 2 quiver gauge theory on x with 
non-zero VEV's for the scalar fields and zero i?-symmetry chemical potentials. We carry 
out this computation in the special case where the background fields are assumed to be 
"commuting" in a sense that conforms to the quiver structure. Furthermore the background 
fields will be taken to be static and spatially homogeneous in order to preserve the >S'0(4) 
isometry of the spatial manifold. The method employed for computing the effective 
potential will be the standard background field formalism. That is, we expand the quantum 
fields about classical background fields and path integrate over the fluctuations, discarding 
terms of cubic or higher order in the fluctuations. The resulting fluctuation operators turn 
out to have a particular tridiagonal structure in their quiver indices. By exploiting the 
vacuum structure of the theory we find that the determinants factorize, leading to an 
expression for the quantum effective action of = 2 U{N)^ quiver gauge theory that 
explicitly displays the Zm structure of the theory. Finally we generalize our results to a 
specific class of field theories that can be obtained as Zm projections of A/" = 4 SYM theory, 
of which J\f = 2 quiver gauge theory is a special case. 

In Section 6 we find the minima of the matrix model of Section 5 in the large N limit in 
a coarse grained approximation. We consider the joint eigenvalue distribution of the scalar 
VEV's and the Polyakov loop and find that the topology of the eigenvalue distribution is 
tied to the Hagedorn phase transition. Below the Hagedorn temperature the eigenvalues 
of the scalar VEV's are distributed uniformly over an S^/Zm and the eigenvalues of the 
Polyakov loop are distributed uniformly over an S^. Thus, the joint eigenvalue distribution 

^We are using a somewhat sloppy terminology here: by 'Polyakov loop' we really mean the holonomy 
matrix of a closed curve winding about the thermal circle and not just its trace. Throughout this paper we 
will use the word to describe both and leave the precise meaning to be determined from the context. 
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is an /Ijm fibered trivially over . We interpret this /Ijm as the emergence of the 
/"Lm factor of the holographically dual AdS^ x /I^m geometry. Above the Hagedorn 
temperature the eigenvalue distribution of the Polyakov loop becomes gapped and is thus 
an interval. The scalar VEV's are now distributed uniformly over an /'Lm fibered over 
this interval, with the radius of the /T^m at any point in the interval proportional to 
the density of Polyakov loop eigenvalues at that point (for fixed TR). The /I^m thus 
shrinks to zero radius at the endpoints of the interval: the topology of the joint eigenvalue 
distribution is an where the Zm is understood to act on the transverse to an 

diameter. Finally we generalize our results to the Zm orbifold field theories discussed at 
the end of Section 5. In particular we find that the geometry of the dual AdS spacetime is 
mirrored in the structure of the quantum effective action in a precise way within this class 
of orbifold field theories. 

In Section 7 we discuss the results we have obtained in this paper and suggest directions 
for future study. In Appendix A further details about N = 2 U{N)^ quiver gauge theory 
are given, some of which the authors of this paper have not found elsewhere in the literature. 
In particular, we write the full Lagrangian density in terms of SU{2)r x U{l)ji invariants. 
In Appendix B we give further technical details of the computation of the quantum effective 
action obtained in Section 5. 

2. M = 2 quiver gauge theory with i?-symmetry chemical potentials 

In this section we review N = 2 U{N)^^ quiver gauge theories on x S'^ with i?-symmetry 
chemical potentials. An introductory review of = 2 quiver gauge theories on 5"^ x 5'^ 
is given in Section 2.1. Details, some of which the authors have not found elsewhere in 
the literature, are deferred to Appendix A. In Section 2.2 we then write up the complete 
Lagrangian density including i?-symmetry chemical potentials. 

2.1 Review of = 2 quiver gauge theory 

N = 2 quiver gauge theory with gauge group U{N)^'^ arises as the world- volume theory 
of open strings ending on a stack of N D3-branes placed on the orbifold C'^/'Lm- The 
gauge theory is thus superconformal |^] with 16 supercharges. It can be obtained as a Z^f 
projection of AA = 4 U{NM) SYM theory as explained in detail in Appendix A. The result- 
ing gauge group is U{N)^ where all the U{N) factors of the gauge group have the same 
gauge coupling constant ^ym associated with them. Letting i = 1, . . . ,M and identifying 
i ~ i -|- M, the field content can be summarized as follows. There are M vector multiplets^ 
{Afj_i, ^i,ip^^i, ipi) where A^^i is the gauge field, ipi is the gaugino, $j is a complex scalar field, 
and ip^^i is the superpartner of We take ^p^ and ip^^i to be 2-component Weyl spinors. 
Furthermore there are M hypermultiplets j, XA,i) X-B,i) where and 

are complex scalar fields and XA,i and XB,i are their respective superpartners which 
we will take as 2-component Weyl spinors. The fields in the z'th vector multiplet all trans- 
form in the adjoint representation of the i'th U{N) factor of the gauge group. The fields 

'^We will use an A/" = 1 notation throughout since this proves convenient. 
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in the i'th hypermultiplet transform in a bifundamental representation of the i'th and 
(i + l)'th factors. More specifically, letting N, denote the fundamental representation of the 
i'th U{N) factor and Nj the corresponding antifundamental representation, (j+i) and 
its superpartner XA,i transform in the Nj (g) Nj+i representation, whereas -B(j_|_i) j and its 
superpartner XB,i transform in the Nj Nj+i representation. 

The field content is conveniently summarized in the quiver diagram in Figure 1. The 
diagram consists of M nodes, labelled hy i = 1, . . . , M with the identification i ~ z + M. 
The i'th node represents the i'th U{N) gauge group factor. Fields belonging to the i'th 
vector multiplet are drawn as arrows that start and end on the i'th node. For the i'th 
hypermultiplet, the fields transforming in the Nj(giNj_|_i representation are drawn as arrows 
that start at the i'th node and end at the (i + l)'th node; the fields transforming in the 
Nj (g) Nj_|_i are depicted as arrows going from the (i + l)'th to the i'th node. 




Figure 1: Quiver diagram summarizing the field content of N — 2 U{N)'^^ quiver gauge theory. 
Each of the black dots (called nodes) represents a U{N) gauge group factor. The nodes are labelled 
by i = 1, . . . , M with the identification i ~i + M . Arrows go from fundamental to antifundamental 
representations of the corresponding gauge group factors. The scalar fields ^j_(i_|_i), and <i>j 

are shown in the figure, whereas the gauge fields and all the superpartners have been left implicit. 



The holographic dual of M = 2 quiver gauge theory was found in |25| to be Type 
IIB string theory on AdS^ x S^/Zm- The quotient S^/Zm is obtained by embedding 
in where the action of Zm is as defined in (|A.1|). The AdS^ space has a radius given 



by R^^g = \/4:TTgs{a'yNM where gs is the Type IIB string coupling. There are also NM 
units of 5-form RR-flux through the AdS^. Due to the orbifold action the volume of the 
quotient /Zm equals the volume of the covering space divided by a factor M where 
the 5^ has the same radius as AdS^. Similarly, there are N units of 5-form RR-flux through 
the /Zm factor which originate from NM units of flux in the covering space. Finally, 
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we note that the Yang-Mills coupling for each U (N) gauge group factor qym is related to 
the Type IIB coupling by ^'yj^ = AngsM. This means that the 't Hooft coupling relevant 
for each factor is A = gyu^ ~ '^'^ds^M. This is the same as the 't Hooft coupling on the 
original NM D3-branes before orbifolding, for which the Yang-Mills coupling was equal to 
4:TTgs. In the following we will often denote the Yang-Mills coupling simply by g. 

The action of = 2 U{N)^ quiver gauge theory defined on x is given as follows. 
To fix our conventions, we set = d^A^, — di,A^ + ig[A^, Ay] and -D^ = d^ + ig[Afj_, •]. We 
will denote the circumference of the thermal circle with /? and the radius of the spatial 
5^ with R. The Euclidean action of = 2 quiver gauge theory on 5^ x S'^ is then 



S 



L 



gauge 



scalar 



~l~ -^ferm) 



(2.1) 



where the gauge boson, scalar field and spinor field Lagrangian densities are given by, 
respectively^ 



C 

-^gauge 
•^scalar 



1 



4 
Tr 



■TrF,,„F, 



D^A D^A + Df,B D^B + D^^ D^^ 



+ R-'^[AA + BB + ^^) + ^g'[[A,A] + [B,B] + [$,$] 



(2.2) 



2gH\[A,B]Y + \[AM +\[BM\ 



(2.3) 



Aerm = « Tr [XA T^,Df,XA + XB T^D^^XB + V' Tf^D^^^J + V'* T^D^^i^^ 

-^{[A,xb] - [B,xa\) - 4^^{[A,xa\ + [B,x^]) 
+ Xa(P,V^] - [B,i>]) + xb{[A,^] + [S,^]) 

- tl;{[A,XB] - [b,xa]) - i^ilAxA] + [B,xb]) 

+ xa[^,xb] - Xb[^,Xa] + - V'*!^,?] 

+ xa[^,xb] - xb['^,xa] + ^^1^,^^] 



(2.4) 



The traces are taken over the NM x NM matrices. The spinor fields xa,Xb,''P<^,'>P are 
undotted 2-component Weyl spinors. We define = (l,icr). The operator D/^ is defined 

by ipiD^4'2 = \{''PiD^ip2 — {D 11)2) ■ It is implied that the fields A^B,(^^A^ etc. take 
the orbifold projection invariant forms given in Eqs. ( [A.15D -( pV.16| ) and ( |A.31| )-( [A.32 ). Note 
that the scalar fields are conformally coupled to the curvature of the spatial manifold 5^ 
through the term R~^ Tr (^AA + BB + in (|2.3| ). This effectively induces a mass for the 
scalar fields. 



*Note that for all fields, including the Weyl spinors xa,XB)^*)^i the bars denote the Hermitian con- 
jugate, not the complex or Weyl conjugate. E.g., (xA)a/3 = (x^)Ja where a, {3 are gauge group indices and 
the * denotes complex conjugation. Furthermore, in the third line of Eq. (2.3), the notation means, e.g., 
\[A,B]\'' = [A,B]\A,B]. 
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The orbifolding breaks the i?-symmetry group SU{4) of = 4 SYM theory into 
SU{2)ji X U{l)ji. As described in Appendix A, $ is associated with the zi direction of 
which is inert under the action of the orbifold group Zm, while A and B are associated with 
Z2 and zs respectively. The U{l)ji factor corresponds to the transformation zi — > e'^zi and 
therefore acts on the $ fields by multiplying phase rotations. The A and B fields have zero 
charge under U{1)r. The SU{2)ji symmetry acts on the A and B fields and their Hermitian 
conjugates. In fact, {A, B) and {—B,A) form SU{2)fi doublets. Furthermore {ip,ip^) and 
{—tp<^,ip) are SU{2)ji doublets whereas xa and xb have zero charge under SU{2)ji. The 
gauge field is not charged under SU{2)ji x U{l)ji. We summarize the i?-charges in Table 
A. 

2.2 Lagrangian density with i?-symmetry chemical potentials 

Given any non-Abelian symmetry group G, one can introduce chemical potentials conjugate 
to the generators of a maximal torus of G. In this section we will consider the case where 
G is the i?-symmetry group SU{2)ji x U{1)r oi M = 2 quiver gauge theory. The maximal 
torus is U{1) X U{1). We will denote the Cartan generators of U{1)r and SU{2)ji by Qi 
and Q2, respectively, and the corresponding chemical potentials by /xi and fi2- For the 
U{1) factor of the maximal torus that corresponds to U{1)r the eigenvalues of the Cartan 
generators can directly be read off from Table A. For the U{1) C SU{2)r we choose as 
a basis for the Cartan subalgebra the diagonal generator cr^ so that the SU{2)ji doublets 
will have well-defined charges under U{1). (We choose az rather than ^cr^ as the generator 
Q2 because we require e'*^^^ to be invariant under 9—^9 + 2ir. Setting Q2 = (Tz we have 
^iQ2e — diag(e*^, e~*^) which is clearly invariant.) Therefore the charges under the maximal 
torus C/(l) of SU{2)p> will be 2 times the SU{2)ii charges. 
Thus for the bosonic fields, 

(flaQa) Ai,{i+1) = (2-5) 
(/^a<3a) = /'''2 (2.6) 
{^^aQa) = (2.7) 

Ma) A^i = 0, (2.8) 

and for the fermionic fields, 

{lJ'aQa)XA,i = -^IJ'lXA,i (2.9) 

{lJ'aQa)XB,i = -\lJ'lXB,i (2.10) 

{iJ-aQa) i^i = - IJ-2) i^i (2.11) 

(MaQa) = (-5/^1 - 1^2) ■ (2.12) 

The corresponding expressions for the Hermitian conjugate fields are obtained by simply 
changing the signs of the chemical potentials. 

To obtain the Lagrangian density of A/" = 2 quiver gauge theory with chemical potentials 
Ha for the SU {2)r x U {1)r Cartan generators, one makes the following substitution in the 
Lagrangian density 

D^- IXaQa&I^O. (2.13) 
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Below we have written the Lagrangian densities for the fundamental scalar and spinor 
fields of = 2 quiver gauge theory. This will be important for the analysis in the following 
sections in order to distinguish the adjoint from the bifundamental structures. 

The Lagrangian density for the scalar fields with i2-symmetry chemical potentials is 



M 



^scalar " X] ^ f (j+i) + Z^A^j^dj (j+i) - ^^(j+l) - /i25/x0^i,(i+ 



1) 



1=1 



+ Tr 



+ Tr 



X {9fiAi,{i+l) + '^9^fj.{i+l)A,{i+l) - + ^2'^Aj0^i,(j+l) 



+ - A{i-i),i^(i-i),i 



2^2 Tr 



2^2 TV 



2g2Tr 



A,{i+l)^ii+l),i ~ 



B(i+l),i^i - ^i+lB(^i^i) i j ( ^i+i - 



(2.14) 



Here the traces are always taken over the gauge indices of the Nx N matrices. Observe that 
the chemical potentials and ^2 act like negative mass squares for and (j+i), -B(j_|_i) j. 
On a compact spatial manifold such as S^, these terms are balanced by the positive mass 
square terms induced by the conformal coupling to curvature. We immediately observe from 



( 2.14 ) that M = 2 quiver gauge theory on 5^ x 5'^ is well-defined as long as fii, fj,2 < R^^. If 
the chemical potentials exceed this bound, the theory develops tachyonic modes and there 
exists no stable ground state. 

The Lagrangian density for the spinor fields with i?-symmetry chemical potentials is 



M 



i=l 



Tr (^{d^XAi + ^g^^,{^+l)XA~i - ^OXAi^i - IfJ^i^fiOXAi) XA,i 
+^ Tr (xBjT^{d^,XB,i + ig\{i+i)XB,i - igXB,iA^i + ^fJ'iS^,oXB,i 

Tr Udf^XB^i + igAf^iXB^i - WXB^iA^ii+i) - h^^l^^^oXB^) ^"^ XB-. 
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Tr {^{d^^i + iglA^i^ilji] + (i/ii - /^a) (5;,oV^ "^/^ 
+^Tr (^V<i>,i Tf^{^f^^p^^i + ig[Af,i,ipq,^i] - (i/ii + /^a) (5/,oV'*, 

+ e'^^lxAK^Oc + ^'"^{XAU {Xi)c{Xi)d} 

+ e''^{xi?,i(Ai)c, (WT)4 + e''^{xB,h (jG)ci\)d} 

- e'"^{xB7tQ^i)c, {Xi)d} - e^'^{xB7i, (Xi+i)c(\)d} 

+ {xA,i^i+i, XB,i} + {xA;7^j, 

- {xB,i'^i, XA,i} - {xBTiWn, xaa}) I • (2.15) 

Here the traces are always taken over the gauge indices of the NxN matrices. Note that the 
potential part of the Lagrangian density has been written in terms of the SU{2)pi doublets 
given in Eqs. ( |A.43| )-( A.44 ) for notational simplicity. 



Finally, as the gauge fields have zero charge under SU{2)r x U{l)ii, the gauge field part 
of the Lagrangian density is unaffected by introducing the i?-symmetry chemical potentials. 
Nonetheless, we give the result here for convenience: 

I 

-^gauge = ^ '^''^ Pfiu^fiu (2-16) 
i=l 

where of course = — d^A^^ + ig[A^^,A1] and the trace is taken over the gauge 

indices of the NxN matrices. 



3. Zero-coupling limit and the matrix model 

The matrix model we will consider is defined by integrating out the fluctuations of the 
quantum fields. In Section 3.1 we therefore first give a brief description of how to compute 
the one-loop quantum effective action with non-zero chemical potentials conjugate to the 
i?-charges. The details of this computation are well-described in the literature (see, e.g.. 
Appendix A of [||). In Section 3.2 we then proceed to construct the matrix model out of 
the 1-loop quantum effective action. 

3.1 One-loop quantum effective action 

The partition function for the grand canonical ensemble has the path integral representation 
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with >Cgauge ! >Cscaiar and >Cfemi being the Lagrangian densities with i?-symmetry chemical 



potentials given by Eqs. ( 2.16|) , ( p. 14 ) and ( 2.15 ), respectively, and where the measures 



T>Afj_, V(j) and Vip are the products of the measures over all the gauge fields, scalar fields 
and spinor fields, respectively. We will obtain an effective action from this expression by 
taking the free limit 5 ^ of the tree-level action. However, since the theory is defined on 
a compact spatial one must impose the Gauss law constraint that all states be gauge 
invariant. We perform the projection onto gauge invariant states by using Aqi as a Lagrange 
multiplier, 

^i^iix) — > A^i{x) + 6f,oai/g (3.2) 

where A^i integrates to zero over 5^ x and Oj are constant Hermitian matrices which by 
gauge invariance can be assumed diagonal, Oj = diag{ql, . . . , qf). To obtain the correct zero 
coupling limit one inserts the decomposition (3.2) into the action given through ( 2.14| )-( [2?T^ ) 
and then takes the 5 — > limit. 

As the quantum fields are defined on x one decomposes them into Fourier modes 
on and 5^ spherical harmonics. More specifically, let r denote the direction along the 
. We will use the convention that any field (/> defined on x has the Fourier mode 
decomposition 

00 



{t,x) 



fc=— 00 



(3.3) 



where the quantized Matsubara frequencies are ujk = for bosons and ujk = ^^^^^^^ 
for fermions giving, respectively, periodic and antiperiodic boundary conditions around 
the thermal circle.^ One then decomposes the spatial components of the gauge field into 
spherical harmonics on by writing them as a sum of a transverse (i.e. divergenceless) 
vector field and a longitudinal vector field VFi where Fi is a scalar function. That is, 
for /c = 1, 2, 3 we decompose 

Af = (Ai)'' + (VFi)^ (3.4) 
and insert the expression on the right hand side into the action given through (|2.14| )- (|2.1(^ ). 



quantum field 


eigenvalue 


notation in text 


degeneracy (Dh) 


transverse vector A"*- 


-{h + l)'^R-^ 




2h{h + 2) 


longitudinal vector VF 


-h{h + 2)R~'^ 




{h + lf 


real scalar A^, (j) 


-h{h + 2)R-'^ 


-A^ 


{h + lf 


Weyl spinor ip 






h{h + l) 



Table 1. Eigenvalues and corresponding degeneracies of the spatial Laplacian = + + 
for various quantum fields defined on . Here R denotes the radius of . The irreducible represen- 
tations of the S'0(4) isometry group are labelled by the angular momentum h which has the range 
/i = 0, 1, 2, . . . for all the fields except for the longitudinal vector field where h starts from 1. 

^However, for the Fadeev-Popov ghosts the boundary conditions are taken periodic. 
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The quantum effective action F = — In Z is defined by integrating out all fluctuating fields 
(cf. (|3.1|) ), leaving an expression that only depends on the zero mode Oj. It is convenient 
to express T as a functional of the holonomy matrix of a closed curve wound around the 
thermal circle, i.e. Ui = e*^°' after decomposing the gauge field according to ( |3.2| ) and tak- 
ing g ^ 0. By performing the traces over the Matsubara frequencies and over the angular 
momenta /i, with appropriate eigenvalues of the Laplacian on and the associated de- 
generacies (cf. Table 1) one finds the following expression for the quantum effective action 
in terms of the variables x = and yj = e^'^J 



nui] 



M oo 

EE 

i=l 1=1 



M oo 

EE 

1=1 1=1 




[y{ + yi') 



1/2 . -l/2\ 

Vi +yi ) 



X^+X^^ \ , I (-1) 

^2 + ^2) + 



l+l 



X (Trt/iTr[/^\ + TrC/r'Tr[/i+i 



(2/2 + ^2 



{l-x^Y 



Tr Ul Tr f/r' 



[vi +yi 



2 -1/2-^ 



(3.5) 



Note that the adjoint holonomy factors come from the vector multiplets (j4^i, ^i,tp^,i,tpi) , 
and the bifundamental factors come from the hypermultiplets -B(j_|_i) j, XA,ii XB,i) ■ 

For later convenience we define here the total single-particle partition functions for the 
bosonic and fermionic sectors of the vector and hypermultiplets: 



z^d[x;yi,y2) 

z^dix;yi,y2) = 

Zhi{x;yi,y2) = 

zu{x;yi,y2) = 



6x 



(1-X^3 



Oq^'J ,71 _|_ rv''^ 

+ (yi + 2/1 ') 



)3 (l-x)S 



(1-x) 

X + 

(1-X)3 
23.3/2 

(1-X)3 



/ 1/2 , -l/2w , -1^ 

^{yi +yi ){y2 + y2 ] 



(2/2 + ^2 ^) 

/ 1/2 , -l/2\ 

[yi +2/1 ) ■ 



(3.6) 
(3.7) 
(3.8) 
(3.9) 



These results are consistent with Ref. Eqs. (3.17)-(3.18), where the summation over 
representations is taken to run over the adjoint and the bifundamental representations, and 
the charges Q are taken as /3 times the Cartan charges Qi,Q2 given implicitly through 

dD-dUl). 

3.2 The matrix model 

The matrix model we will consider is defined by the partition function 

M 



Zmm = J fl [DUi\ exp {-T[Ui]) 
1=1 



(3.10) 
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where T[Ui] is given in (|3.5| ). It is convenient for taking the continuum Hmit to rewrite T[Ui] 
directly in terms of the zero modes a,. To simphfy the notation, define the rescaled zero 
mode ai = f3ai so that Ui = e*"» . Hence 

Zmm = I II [Dai] exp(-Y, (Kd(ar " + Vhi{aT - c^+i))\ (3.11) 

i=l \ m^n / 

where the adjoint and bifundamental potentials are, respectively 



V^{e) = - In 



sm 



A oo ^ 



oo ^ 

= ln2 + (l - zfd(^';2/U2) - {-l)'^'zl^{x'-J,J,)) cos(W) (3.12) 
1=1 

CO 

Vu{0) = - Ey {^u{x';y{,yi) + {-iy^'z^,{x';y[,yi)) cos(/0) . (3.13) 
1=1 

We will now take the continuum limit — > cxd. It is convenient to introduce eigenvalue 
distributions pi{9i) proportional to the density of the eigenvalues e*^' of Ui at the angle 
6i € [— 7r,7r]. Here pi must be everywhere non-negative, and we choose its normalization so 
that for any fixed i 

dOi pi{9^) = 1 . (3.14) 



Furthermore we define the Fourier modes of pi and T4d and Vu- 

p\ = r d9iPi{9i)cos{W^), Kd= r d9V,a{9)cos{l9), V^,= F d9 Vu{9) cos{W) 

J —7T J —TV J —TT 

(3.15) 

SO that, assuming pi^ V^d? to be even functions, we have the Fourier expansions 

^ oo ^ CXD ^ OO 

P^iO = -Y.A ^-d(C) = - E ^ad cos(ZC), Vu{C) = - E ^bi cos(ZC). (3.16) 

«=o (=0 «=o 

The continuum limit is obtained by making the substitution^'' 

^[...] ^ / (3.17) 



n=l 



Furthermore we replace the path integral measure \Vai\ — > \T>Xi\ ■ Thus, in the continuum 
limit the path integral of the matrix model takes the form 

„ M 

/ n [^^^] (-'^mmW) (3.18) 



i=l 

where the action for the eigenvalue distribution functions p is 

^^t2 M oo 

vr 



Smm[p\ = —^^{{p'ifvUT;pi,P2) + pipi+X-AT; pi,P2)) ■ (3.19) 

i=l 1=1 

To summarize, the matrix model under study is defined by Eqs. ( |3.18| )-( ]3!T9| ). 



"'Here it is implied that the content of the brackets \ ■ ■ ■] carries an i label. 
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4. Phase structure 



The term — In |sin (|) | in the adjoint potential ( p.l2 ) originating from the change of measure 



is a temperature-independent repulsive potential. On the other hand, the remaining parts of 
the adjoint and bifundamental potentials ( ^.12|) -( 3T^ ) provide an attractive force^^ which 



grows from zero to infinite strength as the temperature is raised from zero to infinity. One 
would therefore expect that at low temperatures, the stable saddle points of the matrix 
model are characterized by the eigenvalues of the holonomy matrices Ui spreading out 
uniformly over the unit circle, whereas at high temperatures the attractive potential causes 
them to localize 

4.1 Low-temperature solution and phase transition 

We now consider the saddle points of the matrix model action ( p. 19 ), 



= ^ = ^{2p\vU + {pU + Pl.)y^) ■ (4.1) 
For M > 2, this condition translates into M linear equations in M unknowns: 

2piKd + {Pli+Pi+M = 0. (4.2) 

The determinant of this system of equations is generically non-zero, so we find the unique 
solution p\ = 0, corresponding to the flat distribution Pi = Thus we conclude that the 
eigenvalues of the holonomy matrices Ui are distributed uniformly on each of the M unit 
circles. This defines the low-temperature solution of the matrix model. 

The leading 0{N'^) contribution to the free energy computed from the path integral 



( 3.18 ) comes from the action jSmmIp]- However, as = 0, the first non-zero contribution to 
the free energy in this phase comes from a Gaussian integral over the fluctuations about the 
solution /Oj = The free energy is therefore of 0{1) with respect to N, suggesting that 
the theory in this phase describes a non-interacting gas of color singlet states. Furthermore, 
we note that the Polyakov loop W{C) = TrPexp (^ig dx^A^^ has zero expectation value 
since the trace averages to zero in the uniform eigenvalue distribution. In particular, this 
implies that the Zat center symmetry is left unbroken in this phase. Accordingly, we label 
this phase "confining". 

For M > 2 the solution Pi = ^ will be a minimum of the action until we reach values 
of {T;pi,p2) for which 

dpi dp] 

for any fixed I. When the temperature or the chemical potentials are raised above these 
critical values, the flat distribution becomes an unstable saddle point of the matrix model, 
and the model thus enters a new phase which we will discuss in the next section. For now 



= det Hij 



(4.3) 



we note that (4.3) defines a phase transition condition of the matrix model. 



^^The fact that the remaining parts of (3.12)-(3.13) are attractive potentials can be shown following the 
argument in |H|, footnote 32. 
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It will be convenient to express the Hessian matrix in terms of the variables S,i = 2K^j 
and rji = V^-^. Note first that in the special case M = 2, due to the identification i i + M = 



i + 2, the Hessian matrix obtained from ( 3.19 ) takes the form^^ 

H = 




(4.4) 



The determinant factorizes as det-ff = — 4(?7; — ^6) + ICz)- For M > 3 the Hessian 
matrix is a tridiagonal, periodically continued matrix: 



(,1 for j = i 
7]i for j = izizl 



(4.5) 



where, as usual, we make the identifications i ~ i + M and j ^ j + M. The determinant of 
H factorizes as follows^^ 



det 



■■ r] 

ri U 



M 



2vrj 
cos I — ) 



(4.6) 



Thus, the determinant of H vanishes on any of the lines + 2 cos (^^^ rji = for j = 
1, . . . , M. To single out the physically relevant condition for the vanishing of detH we will 
first consider the case M = 12 to gain intuition. For M = 12 the determinant in particular 
factorizes as 



detH 



(4.7) 



where I is fixed. In Figure 2 we have divided the {S,i,r]i) plane into regions where H is 
positive-definite (denoted by +) and where H is indefinite (denoted by — ). 

Thus regions marked by + correspond to a local extremum (minimum) of Smm, and 
regions marked by — correspond to unstable saddle points. In Figure 2 we have furthermore 
marked the region occupied by the J\f = 2 quiver gauge theory matrix model in the low 
temperature phase by plotting (^i,7?i) for (T;/ii,/i2) = (0.1; 0.8, 0.8). For fixed chemical 



potentials, -Zad' -^ad' -^bi' -^bi increase monotonically with the temperature. Therefore, as 
the temperature increases, the dot in Figure 2 will move as indicated and hit the instability 
line "qi = — at the phase transition temperature. 

By the same analysis, for any M > 2 the phase transition occurs at the instability line 
rji = a(M)(^/ where a{M) is the numerically smallest negative slope of the zero lines of the 
Hessian determinant. For all M > 3 we find from (4.6) that a(M) = — ^ (corresponding 
to j = M). For M = 2 we also find a{M) = —\- Indeed, note that for M > 2 the matrix 
obtained by substituting -qi = in Eqs. (4^) and p.5| ) will have a zero eigenvalue (with 

(1, 1, . . . , 1) as an eigenvector) and hence zero determinant. 



^^We omit here, and in the following, the overall factor of 
^^This formula is a special case of (|B.i 



in Eq. (3.1£) for notational simplicity. 
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Figure 2: Regions of positive- definiteness and indefiniteness of H for the case M = 12. Regions 
where H is positive- definite (corresponding to a local minimum of Smm) are marked by +; regions 
where H is indefinite (corresponding to an unstable saddle point of Smm) are marked by —. The 
lines represent the locus o/det H = 0. The physically accessible region of the {^i, r/i) plane is bounded 
from above by the axis and from below by the line of the numerically smallest negative slope. This 
is illustrated by the dot which corresponds to {T; ^1,112) = (0.1; 0.8, 0.8) and I = 1. The arrow shows 
how the dot will move as the temperature is increased, keeping /ii and ^2 fixed. 



The large M limit 

As a consistency check, we can derive that limM^ooCi{M) = — ^ by a different route. We 
take the continuum limit M ^ 00 in the quiver direction. The quiver label i thus becomes a 
continuous angular parameter '& which we take to be 27r-periodic; i.e., we identify ~ i}+27r. 
Accordingly we make the substitutions 



where ' denotes The matrix model action ( 3.19| ) thus becomes 
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Sum [p] 



27r 



(I'd 



{ii + 2r^i){p\^)f-^i{p\^)y 



(4.8) 
(4.9) 



(4.10) 



The Euler-Lagrange equations obtained from this action are those of a harmonic oscillator, 

r]ip\^) + {il + 2r]i)p\^)={) (4.11) 

where / = 1, 2, . . .. Note here that it is the bifundamental contribution in ( 3.19| ) that gives 
rise to the derivative term in ( [1.10| ) and in turn to the mass term for the harmonic oscillator. 



'*The extra prefactor |^ comes from changing the counting measure over i to the measure di). 
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Thus, the harmonic oscillator EOM's in the large M limit is a pure 'quiver phenomenon'. 
Solutions to these equations will become unstable when the tension r = + 2rii) goes 
from r > to r < 0. Thus, for large M, the phase transition will occur when r]i = 
consistent with what we found above. 

We now return to the phase transition condition rji — ol{M)^i. Since -^^j ^i]^; 
are all monotonically increasing as functions of x and < x < 1, the / = 1 condition is the 
strongest. Therefore, the phase transition condition for M > 2 is 

for M > 2 : [z^{x;yi,y2) + z^(x;yi,y2)) + 2{z^^{x;yi,y2) + z^^{x;yi,y2)) = 1 . (4.12) 

Finally, in the special case M = 1 we immediately obtain V!^^ + = from ( p] ) due to 
the identification i~i + M = i + l. Putting 1=1, this is precisely the phase transition 
condition ( 4.12| ). We thus conclude that for any M the phase transition condition is 



{z^{x\yi,y2) + z^{x;yi,y2)) + 2{z^i{x;yi,y2) + z^^{x;yi,y2)) = 1 . (4.13) 

In Figure 3 below we have plotted the curves in the (T, fi) plane obtained from this condition 
for the cases {^1,^2) = (^,0) ; {^1,^2) = (0,^) and {^1,^2) = (M;/^)- For each of these 
cases, the relevant curve defines the phase diagram of = 2 quiver gauge theory as a 
function of both temperature and chemical potential. Note that, as discussed in Section 
2.2, if one or both of the chemical potentials are larger than the inverse radius of the spatial 
manifold 5"^, the theory develops tachyonic modes and becomes ill-defined. Therefore the 
line ^ = 1/R defines a boundary of the phase diagram. 

The phase transition condition ( [4. 131 ) defines a phase transition temperature Th{hi, IJ-2) 
as a function of the chemical potentials. We will refer to Th{^ii-, ^2) as the Hagedorn tem- 
perature of = 2 quiver gauge theory. This terminology will be justified in Section 4.2. 
We remark that the Hagedorn temperature at zero chemical potential is 

Th = ^- — ^ ~ 0.37966 (4.14) 

ln(7-4V3) ^ ' 

in units of the inverse radius of the . This is exactly the Hagedorn temperature for 
= 4 SYM theory (cf. |^, 0). The origin of this fact can be traced to the observation in 
1 30 1 that in the large N limit the correlation functions of = 4 U{N) SYM theory equal 
the corresponding correlation functions of the M = 2 quiver gauge theories obtained from 
orbifold projections. Since our computations rely on perturbation theory (namely, taking 
the limit of the action and then performing Gaussian path integrations), and we are 

furthermore taking the N ^ 00 limit, we should expect that the matrix model defined out 
of the quantum effective action will have the same behavior for the M = 2 quiver gauge 
theory as for the AA = 4 SYM theory. 

Furthermore, for small chemical potentials the Hagedorn temperature is given by 

TH{fil,f^2) = ^ + c(m? + + Cufxf + Cufiinl + C22^t + 0(/uf ) (4.15) 
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where the coefficients are 



/3o = -ln(7-4y3), c = , cn = (4.16) 



Ig ' '^11 864 -627+362^3 J 

_ _§o_ f 1810/30-1045 V3/3o-2896v^+5016 \ — lo_ ( 362/3o -209 V3/3o- 1448v^+2508 '\ / . -, 

•^l^ 216 -627+362v^ ) ' '^^2 iqs -627+362^3 7 ' 




Figure 3: Phase diagram of Af = 2 quiver gauge theory. The outermost curve is the transition line 
corresponding to {^i , — (/^! 0) . It has slope in the neighborhood of the point (T, /i) = (0, 1) . The 
inbetween curve corresponds to {^1,^2) — (0,/i), with slope — In 2 near (0, 1). The innermost curve 
corresponds to {^1,^2) = slope — ln4 near (0,1). The phase transition temperature at 

zero chemical potential is common for the three curves and equals T = — ^^i^j^^-^ ~ 0.37966 as in 
theN SYM case. 



4.2 Solution above the Hagedorn temperature 

As the temperature is increased beyond T > Th, the attractive terms in the pairwise 
potential continue to increase in strength, and so the eigenvalues will become increasingly 
localized. The precise distribution can be determined, following |^, by the condition that 
a single additional eigenvalue added on the i'th circle experiences no net force from the 
other eigenvalues on the circles i — 1, i and i + 1: 

= r dC 2K'd(«i - C)P^{C) + r dC V^,{ai - C) (pi-i(C) + Pi+liO) (4.18) 

J —TT J —TV 
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where V^d and Vu are given in ( ^.12 ) and ( 3.13 ), respectively. This provides M equihbrium 
conditions for the lattice action 



M oo 



Si. 



att 



i=l 1=1 



where 



k = z^,ix';y{,yi) + i-iy+h^,{x'-y[,yi) . 
The exact solution of ( 4.ig| ) was found in ||3l[]. It takes the form 



(4.19) 



(4.20) 
(4.21) 



p,(0) = -(sin^(^ 



sm 



2 / e~a. 



1/2 



J2 Qk COS {{k - l/2){e - Oi)) , i = l...M (4.22) 



k=l 



where 

oo 

The support of pi is [oi — 0Q,aj + 6q]. Here one must impose the consistency requirement 



(4.23) 



Pi 



d6pi{9)cos{n9) . 



(4.24) 



For simplicity the following analysis will be carried out only in the truncated case am>i = 
bm>i = which shares the same qualitative behavior with the general case. For n = 1, the 
consistency condition ( 4.24 ) then becomes 

Pi = ^{^^pI + bip}-i + hpl+i) I ^ (siv? (^%) - sin^ ' cos (|) cos(C + ai) 

= cos ai [aip] + hip]_^ + bip]+^) {2si - sj) (4.25) 

where sf = sin^ (^)" analogy with |^], 6q is determined from Q\ = + 2, leading to 
the M equations 

1 = 2sf{aip} + hpl, + feiPi+i) • (4.26) 
By means of (|4.26| ) , one can rewrite (f4.25|) as the set of M coupled equations 

+ [bipl-i + hpj+i - ai cos ai)pj + cos ai{l - hpj^^^ - bip}^^) = . (4.27) 

If one allows some of the Oi to be nonzero, one finds for any fixed pi,p2 a range of tem- 
peratures above Th{pi,P2) where ( 4.27| ) has no solution such that all Sj, given by ( 4.26 ), 



satisfy < < 1. Thus one must have ai 



ttM = 0. 



With these centers of masses of the eigenvalue distributions, requiring that < < 1 
leads to a unique solution of ( [4.27| ). This solution has all pj equal, as well as all Si equal 
and given by 

,2 f ^0 



sm 



1 



1 



1 

ai + 2bi 



(4.28) 
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With the assumption am>i = &m>i = the exact solution ( 4.22| )- (|4.23 ) thus truncates to 



^^'l-sin^(^)cos(^). (4.29) 

It is immediately clear from ( ^.281 ) that for temperatures above the Hagedorn temperature 
one has 9q < vr; i.e., the eigenvalue distribution becomes gapped. In particular we note that 
the phase above the Hagedorn temperature has unbroken quiver translational invariance; 
i.e. Pi = pi+i- The unbroken quiver translational invariance is expected on more general 
grounds due to the perturbative equivalence between = 4 SYM theory and M = 2 quiver 
gauge theory |0, although it should be noted that in Ref. js^ the gauge theories are 
studied on rather than x S'^ }^ 

Free energy slightly above the Hagedorn temperature 

Using the Hagedorn temperature for small chemical potentials given in ( ^.15| )- p!T7| ) we 
can compute the free energy slightly above the Hagedorn temperature in analogy with 
Defining AT = T — TnifJ^i, fJ-2), we find for < AT <^ 1 the perturbative expansion 

AT 





- /?o^^^(^^? + + Oif^t)^ AT3/2 + 0(AT2) . (4.30) 
High-temperature behavior of free energy 

In the T — > oo limit the pairwise attractive potentials grow to infinite strength, so the eigen- 
values of the holonomy matrices Ui localize to extremely small intervals; i.e. the eigenvalue 
distribution functions will become delta functions, pi{Oi) 6{0i). (This is also clear from 



(4.28) since for T ^ oo one has ai, bi ^ oo and thus 0. The normalization condition 



(I3.14D then imphes piiOi) 6{9i).) Therefore p[ 1, and so from ( ^191) and (|U)-(U) we 



find that the free energy in the T ^ oo limit is 



F = -N^M + + 2/.^) -^M- Vol(5^) . (4.31) 

Here we have applied the polylogarithm regularization procedure described in Appendix 
E of Ref. [^] in order to obtain ( 4.31] ).^^ We note that the free energy scales as N^M as 
N ^ oo. This is to be expected from the orbifold projection invariant form of the fields 
( 1A.15| )- (|A.16| ) and (|A.31|) - (|A.32D , given that the free energy scales as A^^ for = 4 U{N) 



SYM theory for high temperatures in the — > oo limit (cf. Eq. (5.62) of 



^^Note that non-perturbative effects could potentially destroy the quiver translational invariance. How- 
ever, the computation carried out in this section is only valid perturbatively, so we would not expect to see 



such effects. See Refs. ||3^, |3^, |34|, ^ ^6|| for work on non-perturbative equivalence between parent/daughter 
gauge theories related by orbifold and orientifold projections. 

^^Note that there is a minus sign missing on the right hand side of (E.4) for the n 7^ 1 case. 
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The fact that the free energies ([4.3Q|) and i \i.31\) are both of 0{N^M) with respect 
to suggests that the gauge theory in the phase above the Hagedorn temperature de- 
scribes a non-interacting plasma of color non-singlet states. Furthermore, from the fact 
that the eigenvalue distribution ( 4.28D -( ^^ ) is gapped we can immediately conclude that 
the Polyakov loop W{C) has non-zero expectation value as the trace does not average to 
zero in this case. In particular, this implies that the Z^v center symmetry is spontaneously 
broken in this phase. Accordingly, we label this phase "deconfined". Thus, we conclude that 
the phase transition defined by Eq. (f4.13 ) is a confinement /deconfinement phase transition. 
Since furthermore the derivative of the free energy with respect to the temperature is dis- 
continuous at the phase transition temperature TH'(/ii, /i2), we conclude that the transition 
is of first order. Furthermore, cf. |0, |6|, we identify it with a Hagedorn phase transition, and 
Tf/(/xi,/i2) is thus the Hagedorn temperature oi M = 2 quiver gauge theory. 



Twisted partition function 



In analogy with |38|, one may study the twisted partition function for the quiver gauge 
theory, taking the boundary conditions for the spinor fields on the 5^ to be periodic rather 
than antiperiodic. In this case the Matsubara frequencies for the spinor fields will be the 
same as for the bosonic fields, and the twisted partition function Z = Tr(— l)'^e~^^ = 
e-r[C/i] jj^g^y obtained directly from (|3.5D by replacing (—1)'"'"^ — > Following |38 



we choose to exhibit the Zm symmetry of the (twisted) partition function by rewriting the 
adjoint and bifundamental holonomy factors in terms of eigenvectors under quiver node 
displacements i — > i -|- 1. Indeed, define for uj = ^'"''■1^^ ^ 



M 



(4.32) 



Under the quiver node displacement Ul — > ^l+i ^^"^ — *■ ^^^^ Jlj, is an 

eigenvector under the displacement with the eigenvalue iv^ . Writing the holonomy factors 
in terms of 0^, the twisted quantum effective action takes the form 



^ M oo r ^ 
M I 

k=l 1=1 



^ M oo 
k=l 1=1 




1/2 , -l/2\ 

Vi +yi ) 



X UJ 



Tr n[ Tr + Tr Tr 



-k 



(4.33) 



It would be interesting to study the phase structure for the twisted partition function. 



4.3 Quantum mechanical sectors 

Since Af = 2 quiver gauge theory is a conformal field theory, we can exploit the state/opera- 
tor correspondence and map the Hamiltonian H to the dilatation operator D. As a conse- 
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quence, the partition function of thermal M = 2 quiver gauge theory in the grand canonical 
ensemble takes the form 



Z{T- f^^,^,2) = Tin (e-^^+'^^»«') . (4.34) 

Here the trace is taken over the entire Hilbert space 7i of gauge invariant operators. For 
weak 't Hooft coupling A ^ 1, the dilatation operator D can be expanded perturbatively^^ 

oo 

D = Z)o + ^A"/2i^„. (4.35) 

n=2 

We let Q denote the total charge with respect to the Cartan generators of SU{2)r x ?7(1)r, 
Q = Qi + Q2, with /i as the associated chemical potential. Taking A = 0, the partition 



function ( 4.34 ) can be rewritten as 

Z{T; fi) = Ttn exp ( - P{Do - Q) - /3(1 - //)Q ) . (4.36) 



Following j^, we now consider the region of small temperature and near-critical chemical 
potential 

T < 1 , (4.37) 

In this region, the Hilbert space of gauge invariant operators of = 2 quiver gauge theory 
truncates to certain subsectors. To show this, first observe that in the region ( 4.37| ), oper- 



ators with Dq > Q appear with an extremely small weight factor in the partition function 



( 4.36 ) since (3^1. On the other hand, for operators with Do = Q, the weight factor 
is non-negligible precisely because 1 — <C 1. Therefore, the partition function ( [4.36 ) is 
dominated by contributions from operators belonging to the subsector 

no = {O en\ {Do-Q)0 = 0} . (4.38) 

We thus conclude that by taking the near-critical limit 

X — > , xy fixed , (4.39) 

the full Hilbert space Ti of gauge-invariant operators effectively truncates to the subsector 
Ti-Q. We will consider three concrete examples of this truncation below, obtained by either 
turning off one of the i?-symmetry chemical potentials, or by putting them equal. As we 
remark below, the resulting subsectors are in a certain sense quantum mechanical. 

Case 1: The 1/2 BPS sector 

We take (/xi,/i2) = (/^>0), and thus the total Cartan charge is Q = Qi. Taking the near- 
critical limit ( [4.39D of the partition function ( ^.lOD then yields 

M / M 



Z{x;y) 



[ n [^^] E E ^ Tr Ui Tr Ur^ . (4.40) 

i=l \j=l 1=1 ) 



^^This was shown for = 4 U(N) SYM theory in |||, 0. 

^^Recall that in Section 2.2 we defined the generator of the Cartan subalgebra of SU{2)r to be iTz rather 



than ^(Tz so that we have the associated charges Qi,Q2 implicitly given through Eqs. (^.5|)-( 2.12). It is 
these charges we are referring to here, rather than the ii-charges given in Table A. 
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Since the scalar field <I>j has Dq = Q = 1, we therefore conclude that the Hilbert space of 
gauge invariant operators truncates to the 1/2 BPS sector spanned by multi-trace operators 
of the form 

Tr(cI>/^^)TV(cI>^=^)...Tr(a>J). (4.41) 

It is clear that in the near-critical limit ([4.39| ) all operators with covariant derivatives 
decouple. Thus all modes originating from defining a field theory on the spatial manifold 
are removed, and the locality of the field theory is lost. In this sense the resulting subsector 
of the field theory is quantum mechanical. 

Case 2: The SU{2) sector 

We take (/xi,/i2) = (0,/i), and thus the total Cartan charge is Q = Q2- Taking the near- 
critical limit ( [1.39D of the partition function ( |3.10 ) then yields 



Z{x; y) ^ jU [VU.] exp ^ Tr f/i TV . (4.42) 



4 = 1 \j = l 1 = 1 



Since the scalar fields (j+i) and j both have Dq = Q = 1, we therefore conclude 

that the Hilbert space of gauge invariant operators truncates to the SU (2) sector spanned 
by multi-trace operators of the form 



ilTr{zfizE---Z^ll) (4.43) 

i=i 

(i) 

where any letter is one of the scalars (j+i) or -B(j_|_i) j. The subscripts ' — denote 

that the quiver labels on the fields in question must trace out a closed loop on the quiver 
diagram in Figure 1 so as to ensure gauge invariance. I.e., an example of a gauge invariant 
single-trace operator is Tr (vli,{i+i)^(i+i),(i+2)-B(i+2),(i+i)^(i+i),i) • 

Case 3: The SU{2\3)/Zm sector 

We take (/ii,/X2) = (^, /u) and thus the total Cartan charge is Q = Qi + Q2- Taking the 
near-critical limit ( [1.39| ) of the partition function ( p.lOj ) then yields 

'2\ ^ 



» M 

•' i=i 



=1 1=1 \ 

M oo 



+ EE^Trf/iTVf/-\ 

i=l 1=1 



(4.44) 

Since the scalar fields '^'i all have Dq = Q = 1, and the Weyl spinor field 

4>^,i has L>o = Q = |; we therefore conclude that the Hilbert space of gauge invariant 
operators truncates to a subsector spanned by multi-trace operators of the form 



fc 

JJ Tr {w\^Xw!^}^ ■ ■ ■ W^j^J) (4.45) 
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(i) 

where any letter W^ ._^ is either one of the scalars (j^x)i *^«) or the Weyl spinor 
4>^,i- Otherwise, the notation is as explained below ( 4.43|) . 

It would be interesting to study this subsector further and determine its symmetry 
group. This group is presumably a subgroup of the SU{2\2>) symmetry observed in the 
= 4 SYM case |||], and determined by the way the orbifolding breaks the embedding 
of S'J7(2|3) into the full 7^ = 4 superconformal group PSU{2, 2|4). 

In H] the authors considered weakly coupled = 4 U{N) SYM theory on x 5^ 
with i?-symmetry chemical potentials in similar near-critical regions of the phase diagram 
as studied here. It was found that the Hilbert space of gauge invariant operators truncates 
to similar subsectors as identified here, namely the 1/2 BPS sector, the SU{2) subsector or 
the SU{2\3) subsector, depending on which chemical potentials are turned on. Furthermore, 
the analysis in |^ was generalized to small, but non-zero 't Hooft coupling A by utilizing 
the 1-loop correction D2 to the dilatation operator (cf. the perturbative expansion (|4.35| )). 
In the large N limit, D2 restricted to the SU{2) subsector becomes the Hamiltonian of 
an SU{2) spin chain; and restricted to the SU{2\3) subsector it becomes the Hamiltonian 
of an SU(2\3) spin chain. What is remarkable is that in both these cases, the spin chains 



are integrable |41 



|, 40 1, and that the truncated Hilbert spaces can be identified with 
subsectors of the complete dilatation operator of = 4 U{N) SYM theory that are expected 
to be closed to any order in perturbation theory. 

For M = 2 quiver gauge theory, the full dilatation operator along with possible inte- 
grable subsectors is not yet completely settled, so we are not able to immediately generalize 
our results to small, but non-zero 't Hooft coupling A. However, we note that much progress 
has been made in this area. In particular, anomalous dimensions of various operators, the 
anomalous dimension matrix restricted to various subsectors, Bethe ansatze and integrabil- 



ity have been investigated in 1 42, 43, 44, 45, 46, 47, 48, 49, 50, 51, B2, 53 1. 



5. One-loop quantum effective action with scalar VEV's 

In this section we will extend the matrix model for M = 2 quiver gauge theory on x 
in Section 3 to include non-zero VEV's for the scalar fields. To this end we calculate the 
quantum effective action at weak 't Hooft coupling to 1 loop in a slice of the configuration 
space of the background fields. To simplify the calculation we restrict to the case of zero 
i?-symmetry chemical potentials. The potential we compute will be valid within the tem- 
perature range < TR <^ A~^/^. The origin of the bound TR <^ A^^/^ comes from the fact 
that R~^ provides a cutoff on the momentum integrals that appear in the loop diagrams 
that contribute to the effective action. Provided that R~^ is much larger than the inverse 
Debye length, one avoids infrared divergences which would require a resummation of the 
thermal mass of the fields. 

The method employed for computing the effective potential will be the standard back- 
ground field formalism. That is, we expand the quantum fields about classical background 
fields and path integrate over the fiuctuations, discarding terms of cubic or higher order in 
the fluctuations. The background fields will be taken to be static and spatially homogeneous; 
thus, the potential obtained from the computation will be a static effective potential. Fur- 
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thermore, we carry out the computation only in a slice of the configuration space in which 
the background fields are mutually "commuting" in a sense that conforms to the quiver 
structure. 

We now proceed with a more detailed description of the calculation. For convenience 
we first rescale all the fields in the J\f = 2 quiver gauge theory Lagrangian density (as given 
in Eqs. ( [O^ ), ( |AT^ ), (|A1|) and with a factor of ^ym as follows 



-^4>. (5.1) 

We then expand the quantum fields about classical background fields by applying the fol- 
lowing transformations to the Lagrangian density 

— ^ + (5-2) 

^ + (5-3) 

+ (Pi (5.4) 
A^i — > A^i + (5^0 "i ■ (5-5) 

The background fields and Ui are assumed to solve the Euler-Lagrange 

EOM's so that they are the VEV's of the corresponding fluctuating fields. We take the 
background fields to be static and spatially homogeneous, i.e. constant on x S^. This 
is to preserve the SO (4) isometry of S'^ as we will not examine the more exotic phases in 
which the vacuum spontaneously breaks rotational invariance. 

The terms of the Lagrangian density arising after the transformations ( |5.2| )-( |F!5| ) are 
grouped by their order in the fluctuating fields. The terms of zeroth order are grouped 
into a tree-level Lagrangian density. The terms linear in the fluctuating fields combine to 
vanish as the background fields are solutions to the Euler-Lagrange EOM's. We discard 
terms containing fluctuating fields to cubic or higher order. The quantum corrections to 
the tree-level Lagrangian density thus arise from path integrations over the terms quadratic 
in the fluctuations. The result will thus be valid to 1-loop order in the loop expansion. 

It is technically difficult to compute the quantum corrections to the effective potential 
for arbitrary background fields. We will therefore only carry out the computation assuming 
that the background fields satisfy the constraints given below. These constraints are analo- 
gous to requiring that the background fields commute, while at the same time they respect 
the quiver structure of the theory. 

First, the Polyakov loops must "commute" with the scalar VEV's: 

ai+i6(j+i)^j - 6(j+i)^jai = , ai&(j+i),j - ^(i+i),iai+i = (5-6) 

[ai,ct)i] = , ["i,'^ = . 



19^ 



Note that with the redefinition of fields in (5.1), discarding terms of cubic or higher order in the 



fiuctuations is the analog of taking the Qym — > limit in Section 3.1. 
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Second, the scalar VEV's must "commute" among themselves: 



a. 













(5.7) 



Since the zero modes ai,(i+i), (/li and Oi are constant over x S^, the tree-level 

action is obtained from the tree-level Lagrangian density by simply multiplying the volume 
of X S^. After imposing the constraints (|5.6|)-(5.7) the tree-level action takes the form 



^(0) 



Tr ( aj,(i+i) ai,(j+i) + + (pi 4>i 

9ym i-i ^ 



(5i 



We choose an gauge defined by adding the gauge fixing action 



So- 



1 1 



4 



dfiAfj_i + i[ai,Aoi] + «C ( o.(i-i),iA(i-i)^i - Aj,(i+i)aj^(j+i) 



(5.9) 



We will furthermore choose the Feynman gauge ^ = 1 for convenience. The virtue of this 
gauge fixing action is that, using ( |5.6| )-( |5?7| ), it cancels terms appearing in the Lagrangian 
density after the transformations (^^)-([0|) that contain both gauge field and scalar field 
fluctuations. Thus, one can do the path integrations over the gauge field fluctuations and 
over the scalar field fiuctuations separately. 



Specification of the vacuum 

We will restrict to the case where all the zero modes (/"j and ai are taken 

to be diagonal N x N matrices.^'' The most general ansatz satisfying all the constraints 



^"Wheii the VEV's are allowed to be off-diagonal, satisfying the constraints ( |5.6| )-(5.7) along with the 
quiver M-periodicity (i.e., a(i+M),(i+M+i) = o,i,{i+i) etc.) ultimately leads to relations between A'' and M, 
such as iV I M. 
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(|5^)- (15J|) is given by 



diagfe^^J, 



b(i+i),i = diag(e~ 
4>i = 4>i+i 



(5.10) 
(5.11) 
(5.12) 
(5.13) 



If we furthermore require the vacuum to respect the gauge invariance and quiver transla- 
tional invariance of the action along with the quiver M-periodicity, the most general form 
is 



0-1,(1+1) = ^ 
4>i = (t>i+i 



(5.14) 
(5.15) 
(5.16) 
(5.17) 



where w = e^""^'^^ and G Z. This is the vacuum we will adhere to in the computations 
throughout this section. We will find that the expression for the quantum effective action 
is independent of the value of k in ( ^.14D - (^.15|) . 



5.1 Quantum corrections from bosonic fluctuations 

There are radiative corrections to the tree-level potential coming from path integrations 
over the part of the action that is quadratic in the bosonic fluctuations. Below we present 
in a bilinear form the part of the action that is quadratic in the bosonic fluctuations, as 
it appears after being added to the gauge fixing action ( ^.9^ and the Fadeev-Popov ghost 
action, and the constraints ( ^.6P -(^) have been imposed. The path integrals will then be 
Gaussian and can be evaluated easily. 

First we introduce some notation. Define 



Br 



I Jrrin 

M ) mn 
/ {Bl^M)mn 

\ {BM^(M~l))mn , 



2 )mn 



\ (^A/,l)mn , 
/ ($l)mn 

mn 



SO that, e.g.. 



(-^ )mn 



2 jmn 5 



M,l jmn 



and A*^^ 



( (^l,2)n 



\{.Am ,l)nm , 



(5.18) 



(5.19) 



(5.20) 
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Furthermore, we define for fixed m, n the fluctuation operators n™"^ n^", ng" and as 
certain M x M matrices (labelled by i, j = 1, . . . , M) whose detailed form is given in ([B.l|) - 



( |B.4D . Then the part of the action that is quadratic in the bosonic fluctuations (including 



the Fadeev-Popov ghosts Cj,Cj) can be written in the form {k = 1,2,3) 



N 



Sb = ^ I t^^2;vl5l (^(A^'^)^„n™"(A^)„™ + i(9fcF'^)^„n-"(9fcF) 

An -I- (T^^ nmn * 
'2 V '—'g -f^Unm /mn '— 'g ^mn 

+ iA^)mn nr ^mn + (B^)™„ Dg" ^mn + ^mn ) (5-21) 



where, as in Section 3.1, the spatial components of the gauge field have been decomposed 
into a transversal (i.e., divergenceless) part (^^)*^ and a longitudinal part (VFi)^. Thereby 
all the fields have been written in terms of spherical harmonics. The path integrations 
over the bosonic fluctuations (j^.]^), j, and can now readily be done and yield 

the formal expression^^ 



rbosh,ai,(i+i),6(i+i),i,</.i] = ^ TV In det n;;^" + Y l^lndetnr 

m,n=l m,n=l 
N N 

+ Y Trlndetng"+ ^ Tr In det □^'^ . (5.22) 



m,n=l m,n=l 



Here the traces are taken over the Matsubara frequencies and over the spherical harmon- 
ics, and the determinants are taken over the i,j indices of the operators □^'^^ n^"^ n™ 



mn 



and n^". Let us define here for convenience 



n,m — 2 



(((ai,(i+l))nn - ■''(ai,(i+l))mm) )nn " )mm) 

+ ((^{j+l),j)nn - {b(i+l),i)mm) ((^(j+l),i)nn " l^"'' (&(i+l),j)mm) 

nn {(pi)mm) (^{4'i)nn i4'i)mm)^ ■ (5.23) 



Now we apply the determinant formula ( [B.6D to the formal expression ( ^.22 ) for Fbos- Then 
we take the traces over the Matsubara frequencies and over the 5^ spherical harmonics. 



^^We are using a rather sloppy notation here as the term involving is to be interpreted as the 

total contribution from the path integrations over the transversal and longitudinal parts of the spatial 
components of the gauge field, the time component of the gauge field and the Fadeev-Popov ghosts. The 



individual contributions are explicitly written out in (5.24) below 
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labelled by the angular momentum h (see Table 1). This yields the following result 

M Af oo 



1 

2M 



+ 



1 

2M 



i,j=l m,n=l k=—oo 
M N oo 

E E E ^'^>o 1^ 

i,j=l m,n=l fc=— oo 

M N oo 



+ (^-l)ir7E E E Tr.>olnf(^. + (ar-ar'"))' + A^ + ^,,;n,. 



^ i,i=l m,n=l k=—oo 

2 Af TV oo 

+ M E E E ^'h>o In + (a 



nn _ Q,™™ 



Vi 



i j'=l m,n=l fc=— oo 



(5.24) 



Here the first line comes from the path integrations over the transverse part of the spatial 
gauge field, and the second line from the integrations over the longitudinal part. The third 
line comes from integrating over the temporal component of the gauge field and the Fadeev- 
Popov ghosts, contributing with the weights ^ and —1, respectively. Finally, the fourth line 
comes from path integrating over the conformally coupled scalar fluctuations. Note that 
there is an exact cancellation between the contributions of all /i > spherical harmonics 
in the second and third line. As we will see in Section 6, the surviving contribution from 
the h = scalar spherical harmonic will be the dominating radiative correction in the 
low-temperature regime. 

After performing the summations over the Matsubara frequencies and writing out the 
traces over the spherical harmonics with the appropriate eigenvalues of and their 
degeneracies (cf. Table 1) we find 



M N 



rbos = ^E E -/?K.n,^)'^' + 2Er 



i,j=l m,n=l 



1=1 



COS (/?/(a 

Nl/2 



't,r, n,mj 



h=0 



oo ^ 



-/3/((h+l)2iJ-2^,;,,,-.„,^)V2 



cos - 



)) 



1=1 



+ 6^(/i + l)2 /3((/i + l)2i?-2 + 

"^ij; n,m ) 



1/2 



h=0 



oo 

2 ^ £g-/3/((/i+l)"i?-^+^;»„ 



COS [Pl{c 



(5.25) 



where Vij-n,m is defined in (|5.23|) . This is the complete result for the contribution to the 
quantum effective action coming from bosonic fluctuations. 
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5.2 Quantum corrections from fermionic fluctuations 

The fluctuating fermionic fields will also give rise to radiative corrections that can be com- 
puted much along the lines of the bosonic corrections. It is convenient to carry out the 
calculation using = 4 SYM notation for the Weyl spinor fields. The quiver structure 
of the action is taken into account by including appropriate factors ^Ic in the fluctuation 
operator as explained in Appendix B.2. 

The fermionic part of the Lagrangian density can be written in AT = 4 SYM notation 
(cf. ( |A.25| )) in the following bilinear form 

■^ferm = ~2 ^ ^ ^ ^ {i^p)i;mm {^p)i]mn) '^ij I /T~"\"' 1 (5.26) 



where the fluctuation operator D^" is given in (|B.ll ). Taking the determinant of D^"^ as 



explained in Appendix B.2, and taking the traces over the fermionic Matsubara frequencies 

+1 

/3 



Uk = (^^+^)^ QYgj. ^}jg spherical harmonics, one finds the following result 



i,j=lm,n=l h=0 \ 



(5.27) 



where Vij-n,m is defined in (|5.23| ). The factor 4 comes from performing 4 path integrations. 
This is the complete result for the contribution to the quantum effective action coming from 
fermionic fluctuations. 

We conclude that the quantum effective action of AA = 2 quiver gauge theory with 
constant scalar field VEV's satisfying (|5.6| )- (|5.7|) is given by 

Teff = 5^ + Tbos + Tferm (5.28) 

where S^^^ is the tree-level action 

^(0) ^ ^n_pn, ^ ^ (a,,(,+i))„„(^^a^)n„ + (6(.+ l),i)nn (^(«+l),i)nra + (0i)nn(0j) nn I ■ 

(5.29) 

and Fbos and Fferm are given in ( 5.25| ) and ( ^.27| ), respectively, with Vij-n,m given in ( 5.23 ). 



Note that the tree-level potential ( |5.29| ) is attractive, whereas the 1-loop quantum cor- 
rections in ( |5.25| ) and ( ^.27| ) are repulsive. As we will see in Section 6, the competition 
between an attractive and a repulsive part of the potential will cause the equilibrium con- 
figurations of the eigenvalues of the scalar VEV's to be hyper surf aces. 

5.3 Generalization to other Zj\f orbifold fleld theories 

The computations in this section and in Appendix B can immediately be generalized to field 
theories obtained as Zm projections of = 4 [/ {NM) SYM theory where the action of Z^f 



- 30 - 



is that in (AJ) with u! replaced by for p G Z. For these theories^^, the quantum fields 
must satisfy the Zm invariance conditions obtained from ( A.14|) and (|A.3Q ) by replacing 
io ujP. In turn, the fields will take Zjvf projection invariant forms analogous to ( A.15|) - 
( A.l(^ ) and ( |A.3lD -( |A.32 ), except that the bifundamental fields will have non-zero entries 
on the p'th super- or sub-diagonal. That is, A and B will have the non-zero entries (j+p) 
and B(^i^p^^i, respectively, and analogously for the respective superpartners XA and xb- 
As a result, the fluctuation operators n'^n ^ ^mn ^ ^mn ^mn (^).(^) and A^^ in 
( B.18|) will have non-zero entries on the p'th super- and sub-diagonals. Therefore, using the 



generalized determinant formula^^ (where 



det 



Zl Z2 Z3 
ZM Zi Z2 
ZM-1 ZM Zi 



\ Z2 Z3 Z4: ■ 



ZM \ 

ZM-1 
ZM-2 

ZI J 



M 



(zi + L0^Z2 + io'^^Z3 



+ u;^''-'^^zm) 



(5.30) 



we see that the fluctuation determinants factorize as in (|B.q ), with io^^ replacing loK We 
conclude that the quantum effective action of these more general Zm orbifold field theories 
is given by the expression ( ^.28 ) where S^^^ is given in ( 5.29| ) and Tbos and Tferm are given 



in ( 5.25| ) and ( ^.27 ), respectively. The only change is that Vij-n,m now takes the form 



i+p),i)nn — {b(i+p),i)mm) {{b{i+p),i)nn — 
nn isi^ijmm) {(,4^i)nn i,4'i)mm ) • (5-31) 



6. Topology transition and emergent spacetime 



In this section we will find the solutions minimizing the effective potential computed in 
Section 5 (given in ( 5.28 ), ( 5.29| ), (|5.25|) , ( ^.27 ) and (^.23|) ) within the temperature range 
< Tii < A"^/2_ stress that, since the effective action of Section 5 is only valid within 
a sector of constant background fields satisfying ( ^.6D -(^), the minima we find in this 
section are not the absolute minima of the gauge theory, and the phase transitions within 
this sector of background fields do not necessarily extend to phase transitions in the full 



gauge theory (cf. |14|). Nonetheless, we will see that the matrix model of Section 5 exhibits 
some interesting dynamics. 

The resulting distributions of eigenvalues will preserve the SU (2) x U{1) /^-symmetry of 
J\f = 2 quiver gauge theory. As in Ref. |13| we believe that due to the preserved i?-symmetry. 



the minima found here are indeed the global minima of the effective action (within the sector 
of constant "commuting" VEV's). The key observation needed for obtaining the solutions 
is that both in the low-temperature regime and above the Hagedorn temperature Th, the 

^^These theories have also been considered in, e.g., Refs. ps], p^ , po| , liof . 

^^We emphasize that the entries of the M x M matrix in (5.3C) are allowed to be complex numbers. 
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eigenvalue distributions for the scalar VEV's and the Polyakov loop can be solved for 
separately. As we will see, the Hagedorn transition causes a change in the topology of the 
joint eigenvalue distribution when the temperature is raised above Th- 

6.1 Low-temperature eigenvalue distribution 

For temperatures low compared to the inverse radius of the (i.e., TR <C 1), one can con- 
sistently discard terms in the quantum effective potential that are suppressed by Boltzmann 
factors^^, and so one obtains the following low-temperature limit of the effective potential 

rri?<l = 2 ^ ^ ( i'^i,{i+l))nn{aiji+l))nn + (&(j-|-l),i)rin(&(i-|-l),i)nn + {(pi) nni'P'i) nn) 

9ym ■ - - ^ ^ 



i=l n=l 



^ M N 



[2 (^((ai^(j+l))„„ - W (ai,(i+l))mm) ((«j,(i+l))nn - 'i^-' (ai,(i+l))mm) 



2M 

i,j=l m,n=l 



+ ~ {b{i+l),i)mm) ((fe(i+l),j)rm - W (^(i+l),j)mm) 

+ {{(Pi)nn - {4'i)mm) {{4>i)nn " i4>i)mm) ) • (6.1) 



We observe that the eigenvalues of the Polyakov loop are not coupled to the eigenvalues of 
the scalar VEV's. Therefore, for low temperatures, the distribution of the Polyakov loop 
eigenvalues will be the same as in the case with zero scalar VEV's treated in Section 4. 
Thus, we immediately conclude from Section 4.1 that the eigenvalues e*"*"" of the Polyakov 
loop (for i fixed) are uniformly distributed over for any temperature below the Hagedorn 
temperature. Note that for a uniform distribution of the angles a"", the terms multiplied 
by Boltzmann factors in ( 5.25| ) and ( |5.27| ) vanish exactly. Therefore we can consistently 
discard these terms as long as the temperature is below Th- 

In order to find the minima of ( |6.1| ) we make the observation that by making the 
identifications 

- (6.2) 

b{i+i),i = t^^+i),* (6-3) 
(pi = (pi (6.4) 

and applying them recursively to ( |6.ip , the low-temperature effective potential reduces to 

rT/?<l = 2 ( («j,(i+l))n?i(ai,(j+l))nn + (^(j-|-l),j)nn(&(j+l),j)nn + {(pi)nn{(pi)nn 

^ M N 

[2 ^((aj^(j+l))„„ — (aj^(j_l_l))mm) ((aj^(j4.1))n7i — (ai,(i+l))mm) 



2 

2=1 m,n=l 



+ ((^( {b{i+l),i)mm)[{h{i.+l),i)nn - (^(i+l),*)mm) 

+ {{(pi)nn — {(pi) mm) {{(pi) nn {(pi)mm ) ■ (6.5) 



^''We will verify a posteriori that this procedure is valid for all temperatures below the Hagedorn tem- 
perature. 
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It is important to note that the identifications ( |6.2| )- (|6.4|) correspond uniquely to the efi'ec- 
tive potential. That is, if one replaces uj by in (|6.2D - (|6.3D , the potential ( |6l] ) will not 
reduce to (^) for general M. To see this, note that, since all M powers of uj appear in 
( |6.1| ), the order of u)'^ must be M . Thus we must have gcd{q, M) = 1 for all M which implies 
q = l. 

We now proceed with finding the minima of (|6.5|). These will be minima of ( |6.lD where 
the identifications (|6^)-( |6^ ) have been made. It is convenient to introduce the dimension- 
less variables 

{Oi)n = ii{ai)nn , (6.6) 
izi)n,l = Pi4>i)nn, (^j)n,2 = P {o-i, {{+!)) nn , izi)n,3 = /3(&(i+l),j)nn (6.7) 

and 

iZi)n = {{Zi)n,l, {Zi)n,2, (^On.s) (6.8) 

SO that {zi)n G C'^ for fixed i and n. Furthermore we introduce a norm on defined by 

/ 3 \i/2 
\\w-z\\ ^ (6.9) 

where | • | denotes the modulus. Written in this notation, ( |6^ ) takes the form 

^TR<^1 = ^-^^^\\{Zi)nf - Yl U^'^n - {Zi)m\\ ■ (6.10) 

We will now take the continuum limit — > oo and describe the eigenvalues of the 
Polyakov loop and the scalar VEV's by a joint eigenvalue distribution pj(0j, Zi) proportional 
to the density of eigenvalues at the point {9i,Zi) (for some fixed i) and normalized as 
J dOid^Zi pi{6i, Zi) = 1. The continuum limit is obtained by applying the substitution 



1 ^ r 

jyEt---]^ deid^zipiiei,zi)[---] (6.11) 

n=l 



in analogy with ( 3.17| ). Here it is implied that the content of the brackets [• • •] carries an 
i label. In the continuum limit, the equation of motion for Zi obtained from ( 6.10 ) reads 

^z.= [ rfV,p.(z'J^i^. (6.12) 

JDi \\Zi-Z-\\ 

Here pi{-) is defined as the average Pi{zi) = f^^d6ipi{6i, Zi), and Di C denotes the 
support for pi. The solution to ( |6.12p is given by the eigenvalue distribution 

f X ^(ll^ill -rj) .„^„. 

where the radius is given by 

A/3 1024 

y/2TT^R 945 ^ ' ' 
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as can be checked straightforwardly. That is, (|6.12|) is satisfied for any when the eigenval- 
ues are distributed uniformly over an with the radius j^J^ ). Since (|61q|) was obtained 



from the low-temperature effective potential ( |6.1j ) by making the orbifold identifications 



p^)-(|6.4|), we thus conclude that the minimum of (|6.1l) is a uniform distribution of the 



eigenvalues of the scalar VEV's over /Zm where the action of is precisely as in ( A.l ) 



This is consistent with |29], as one should expect in the low temperature limit where thermal 
effects are small. Since we found that the eigenvalues of the Polyakov loop are distributed 
uniformly over an for temperatures below the Hagedorn temperature, we conclude fur- 
thermore that the joint eigenvalue distribution of the scalar VEV's and the Polyakov loop 
is jTLu X in this temperature range. 

It is remarkable that the eigenvalues of the scalar VEV's localize to a hypersurface 
in rather than spreading out over the configuration space. The physical origin of the 
localization is essentially common for the matrix model developed here and the matrix 
model of namely the competition between an attractive part of the quantum effective 
potential, and a repulsive part where the latter is generated by the path integrations. We 
interpret the eigenvalue distribution of the scalar VEV's as the emergence of the jl^u 
factor of the holographically dual AdS^ x jTLu string theory geometry. Finally we note 
that the hypersurface jTLu has the isometry group SU{2) x [/(I), resulting from breaking 
the 5C/(4) isometry via the orbifold identifications (|6^)-( [6^ ). Since this is the full R- 
symmetry group SU{2)ji x C/(1)_r of = 2 quiver gauge theory we believe (cf. ||l3|]) that 
the minimum found here is indeed the global minimum of the effective action of Section 5. 

6.2 Eigenvalue distribution above the Hagedorn temperature 

In the matrix model treated in Sections 3 and 4 where the VEV's of the scalar fields were 
zero we observed that as the temperature is increased above Th ~ 0.38 the Polyakov 
loop eigenvalue distributions open a gap. In this section we will examine how this phase 
transition manifests itself in the general case with non-zero scalar VEV's. 

From the radius ( 6.14|) one in particular finds that for low temperatures ||jZj|| » A, so 



\\z-\\ 

that the tree-level term dominates over the quantum correction by a factor ~ » 1. 
On the other hand, around the Hagedorn temperature Th one finds ~ A, and the 
tree-level term and the quantum corrections come within the same order of magnitude. It 
is therefore natural to re-express the effective potential in terms of the new variables 

(Ci)n,l = ^'^{Zi)n,l , (Ci)n,2 = ^'^{zi)n,2 , {Ci)n;i = ^'^{zi)n,'i ■ (6.15) 

The computations in this section will be valid for temperatures in the range < TR <C 
A~^/^. Since we can no longer neglect the terms multiplied by Boltzmann factors, we have to 
consider the full quantum effective action as computed in Section 5 (given in ( 5.28] ), ( 5.29 ), 



dUI), (U^) and (HI). Once again, we apply the orbifold identifications (|6.2|) - (|6.4D , and 



express the result in terms of the variables ^i,Ci- However, the rescaling with the 't Hooft 
coupling A in ( |6.15D will reorganize the perturbative expansion of the effective potential 
into 

Teff = [6,] + Ar(i) [6,, Ci] + 0{\^) . (6.16) 
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Here the 0-loop term is 



M N oo 



1=1 m,n=l 1=1 



where zf^, zf^, z^, are given in Eqs. Q, Q, 
in this case since we are taking = /U2 = here. 
The 1-loop term in ( |6.16 ) is given by 

n 2 o Ar *^ Af 



cos 



[mi}n - m)m) (6.17) 



Sj), (|3.9D , respectively, and yi = ?/2 = 1 



/3 



j=l n=l 



M N f oo 

-ylE E ||(C)n-(CW|| l + 2j^COs(/ 
1=1 m,?i=l \ Z=l 



(^i)n-/(0i)m) . (6.18) 



From the expansion ( 6.16 ) it is immediately obvious that to leading order in A the Qi are 
unaffected by the Cj- Therefore, to leading order, the eigenvalue distributions of the Qi are 
the same as they were in the case with zero scalar VEV's treated in Section 4. The eigenvalue 
distributions of the scalar VEV's can therefore be found by minimizing r(^)[0j, (^j]. Taking 
the large limit of ( |6.18 ) according to ( 6.11j ) one finds 



Ar2 



/3 



i=l 



"^^E / '^^^«^'C'i'a/'^(^^,c.)/5.(^.,c^)iic-dii. (6.19) 

- (6'i)m) which 



i=l 



Here we have used the identity 1 + 2 YL^i cos {l{Qi)n — l{^i)m) = 27r(^((0i)„ - 
is simply the Fourier expansion of the delta function. 

Now we proceed to minimize the action ( |6.19 ). Since the eigenvalue distributions for 
the Polyakov loop and the scalar VEV's can be solved for separately, the joint eigenvalue 
distribution factorizes: 



P^i0^)6{\\a-n{0,)) 



WCiWHl + {dri/d9ify/^Yol{S^) 
Inserting ( 6.20| ) into the 1-loop term ( |6.19| ) one finds 



(6.20) 



_ 2tt^R 



M „ M „ 

i=i i=i 



i=l 
M 



2^ 



2o2 



C'(3 
47r2i?2 



(6.21) 



- 35 - 



where C = ^945^ • The final term only contributes to the 2-loop order distribution of the 
Polyakov loop eigenvalues and can therefore be ignored. Hence for a minimum we have 

r^{e^) = ^P^{e.) . (6.22) 

As we know from Section 4.2, when the temperature is raised above the Hagedorn tempera- 
ture Th, the Polyakov loop eigenvalue distribution becomes gapped and is thus an interval 
[— ^0) ^o]- The scalar VEV eigenvalues are now distributed uniformly over an /Zim fibered 
over this interval, with the radius of the /"Lm at any point 9i in the interval being pro- 
portional to the density of Polyakov loop eigenvalues at 9i (for fixed TR). The jl^u thus 
shrinks to zero radius at the endpoints ±0o of the interval: the topology of the joint eigen- 
value distribution is an where the Zm is understood to act on the transverse to 
an diameter. Thus, the Hagedorn phase transition manifests itself in the general case 
of non-zero scalar VEV's as a change in the topology of the joint eigenvalue distribution 

In order to understand how the jTLu eigenvalue distribution may be realized in the 
dual AdS spacetime we first need to consider the part of the low-temperature distribution 
5^ X /I^M- The eigenvalues of the Wilson line wound around the thermal circle give the 
positions of D2-branes^^ on the T-dual of the thermal circle in thermal AdS^. As the 
temperature is raised higher and higher beyond Th, the Polyakov loop eigenvalues become 
localized to smaller and smaller intervals. On the AdS side one therefore finds a localized D2- 
brane configuration. It was noted in that a similar localization of D2-branes on a spatial 



circle, at finite temperature, was investigated in [55| where it was observed to produce a 
near-horizon geometry containing a non-contractible . Moreover, it was predicted in |55| 
from supergravity that a 5"^ x 5^ — > S"^ topological transition of a Gregory-Laflamme type 
should take place. In the present case, where the dual spacetime is AdS^ x /I^m, we expect 
the appearance of an ITLm in the near-horizon geometry of the localized configuration of 
D2-branes on the T-dual of the thermal circle. 

We now address the important question regarding the stability of the saddle points 
( 6.13 ) and (|6.20|) against off-diagonal fluctuations. As one may read off from ( ^.21 ), the 

mass of the ij entry of a fluctuating scalar field is \j R~'^ + {{^a)ii — {Va)jj)'^ (for notational 
convenience we here use ifa which are real- valued scalar fields related to the complex scalar 
fields by (|A!9|)- (|A.10|) ). In the saddle point (|6l|), one finds from (^1|) that 

{{Va)ii-{^a)jjf (6.23) 

and so the ratio of masses of an off-diagonal fluctuation to a diagonal fluctuation is 

\/l + 0{X^). For small A this ratio is very close to 1. A priori it thus appears possible for the 
fluctuating fields to have off-diagonal elements in this background, causing the background 
to be unstable. 

Despite this we believe that the sector of constant and 'commuting' VEV's is interesting 
since it may have a connection with the dominant saddle points at strong 't Hooft coupling 



^The D2-branes here are T-dual to the original D3-branes. 
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(where the assumption [ipa-,^b\ = seems natural). It is also worth remarking that at zero 
temperature and strong 't Hooft coupling one finds an distribution of the scalar 

VEV eigenvalues |29|. In a different vein, the topological phase transition j'Lu x — > 



/Ijm provides a natural extension of the Hagedorn/deconfinement phase transition for 
the Polyakov loop eigenvalues studied in Section 4, and it is a tantalizing question whether 
it extends to a phase transition of the full gauge theory at weak 't Hooft coupling. 



Finally it should be noted that the effective action (5.28) has other saddle points in 



which off-diagonal fluctuations can have parametrically large masses. (Indeed, compare with 



| ll[| where these saddle points were associated with the Gregory-Laflamme instability in the 
gravity dual theory.) These saddle points are therefore guaranteed to be stable backgrounds. 
In particular, the effective action (|5.28|) has interesting physics beyond the saddle points 
studied in this section. 

Generalization to other Zm orbifold field theories 

The computations in this section immediately carry over to the more general Z^v/ orbifold 
field theories considered in Section 5.3. In this paragraph we remark on the theory defined 
by letting the action of Z^f be that of (|A.l| ) with uj replaced by u;^ for some fixed p G Z. 
The quantum effective action of the corresponding field theory is obtained from that of 
M = 2 quiver gauge theory by defining Vij-^ n,m to be as given in ( F3l|) . The minima of this 
effective action are found by making the orbifold identifications 

Hii+i) = , - , 4>i = (t>i- (6.24) 

The resulting expression for the effective action is then precisely the same as in the case of 
M = 2 quiver gauge theory treated in this section, and the conclusions carry directly over. 
In particular, having made the orbifold identifications (|6.24 ), one finds the low-temperature 



joint eigenvalue distribution x and the high-temperature distribution . Alterna- 
tively, the joint eigenvalue distributions are jTLu x and /I^m, respectively, where 
the action of Zm is precisely the orbifold action defining the Zm orbifold theory. It is im- 
portant to note that the orbifold identifications ( 6.24D correspond uniquely to the quantum 



effective action of the field theory. Indeed, assume that we make the identifications ( |6.24D 
with some u;"^ replacing uj^ . In order for the quantum effective action to reduce to an ex- 
pression involving norms on we must require uj'^ to have the same order as ujP . That 
is, we must have VM G N : gcd(g, M) = gcd(p, M) which implies q = p. Identifying the 
above S^/Zm distribution with the S^/Zm part of the holographically dual AdS^ x S^/Zm 
spacetime, this shows in particular that, within this class of Zm orbifold field theories, the 
geometry of the dual AdS spacetime is mirrored in the structure of the quantum effective 
action in a precise way. 

7. Discussion and conclusions 

In this paper we have investigated different aspects of the phase structure of AA = 2 [/ (N)^^ 
quiver gauge theories. We have set up a matrix model for M = 2 quiver gauge theories on 
5^ X with chemical potentials conjugate to the i?-charges. We then found the stable saddle 
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points of the model as a function of temperature and chemical potentials. More specifically, 
we identified a low and a high temperature phase separated by a threshold temperature 
Tf/(/x) which marks a Hagedorn/deconfinement phase transition. The condition of stability 
of the low-temperature saddle point was translated into a phase diagram oi N = 2 quiver 
gauge theory as a function of both temperature and chemical potentials. We observed that 
in regions of small temperature and near-critical chemical potential the Hilbert space of 
gauge invariant operators truncates to the SU{2) subsector, or to a larger subsector whose 
symmetry group has yet to be determined. More specifically, we found the SU (2) subsector 
when the chemical potential corresponding to the SU (2)^^ factor of the i?-symmetry group 
SU{2)r X U{1)r is turned on, whereas the larger subsector emerged from turning on both 
chemical potentials and setting them equal. 

We then developed the matrix model oi M = 2 quiver gauge theory in a different 
direction, allowing non-zero VEV's for the scalar fields, but setting the i?-symmetry chem- 
ical potentials to zero. We did this by computing a 1-loop effective potential for constant 
and "commuting" VEV's, valid at weak 't Hooft coupling and in the temperature range 
< TR ^ A^^/^. We furthermore obtained the effective potential for more general Z^/ 
orbifold field theories by an immediate generalization. Then we found the equilibrium con- 
figurations of the eigenvalues of the Polyakov loop and the scalar VEV's. The eigenvalues 
of the scalar VEV's localize to a hypersurface in due to a repulsive part of the effective 
potential of a Vandermonde type, originating from the quantum corrections. We found that 
at the Hagedorn temperature the topology of the joint distribution of the eigenvalues un- 
dergoes a phase transition S"^ x S'^ jTLu — > jTLu- Finally, we identified the /Ijm part 
of the low-temperature eigenvalue distribution as the emergence of the jTLu part of the 
holographically dual geometry AdS^ x /'Lm- It should be noted, though, that the latter 
is a dominant geometrical saddle at strong 't Hooft coupling while the "commuting" saddle 
found from the effective potential at weak 't Hooft coupling is not an absolute minimum 
|14, 13]. Extrapolating this identification to high temperatures, we furthermore note that 
the dual spacetime interpretation of the high-temperature is at present not entirely 

clear. We have also generalized the analysis to a class of Ijm orbifold field theories, thereby 
finding that the geometry of the dual AdS spacetime is similarly mirrored in the structure 
of the quantum effective action in a precise way. 

There are several interesting future directions to pursue. It would be interesting to 
investigate other vacua of = 2 quiver gauge theory which preserve less i?-symmetry. In 
particular, such vacua could prove important when the matrix model with non-zero scalar 
VEV's developed in this paper is extended to include i?-symmetry chemical potentials. 

One could also consider the gravity duals of the phase transitions studied in this paper. 
In particular, the effective potential we computed in Section 5 can be used to study the 
manifestation of Gregory-Laflamme instability from the weakly coupled gauge theory point 
of view 1 11]. This would proceed along the lines of |11] where, above a critical temperature 
Tc 3> Th, the effective potential computed for AA = 4 SYM theory was observed to develop 
new unstable directions along the scalar directions accompanied by new saddle points which 
only preserve an 50(5) subgroup of the global 50(6) isometry group. This phenomenon 
was identified as the weak coupling version of the Gregory-Laflamme localization instability 
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of the small AdS^ black hole in the gravity dual of the strongly coupled gauge theory. 

Furthermore, the results obtained in this paper can be applied to compute the Polyakov- 
Maldacena loop [^] at weak coupling. It can also be computed at strong coupling following 
|57], so this is an interesting object to compare at weak and strong coupling. 

It would be very interesting to study the subsectors of the Hilbert space of gauge invari- 
ant operators that we identified in Section 4.3 in further detail, in particular to determine 
the symmetry group of the subsector corresponding to turning on both chemical potentials 
and setting them equal. A further point of particular interest would be to examine whether 
these subsectors are closed under the action of the full dilatation operator in analogy with 

More generally, these results could prove useful to further investigate the corresponding 
spin chain for the N = 2 quiver gauge theory. 



Another direction to pursue would be to examine, following [EG, 55], whether the Hage 



dorn temperature of = 2 quiver gauge theory in the near-critical limit combined with 
the triple scaling limit of |46| can be matched to the Hagedorn temperature of Type IIB 
string theory on a pp-wave with a compact spacelike circle. This would involve computing 
the spectrum of a certain subsector of the SU (2) sector of gauge invariant operators which 
are dual to strings that wind about and/or have oscillators in the compact direction. This 
might require finding novel Bethe ansatze since the ground states of the spin chain govern- 
ing the truncated SU{2) sector appear to be inherently different from, say, ferromagnetic 
ground states. 
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A. Detailed description of AT = 2 quiver gauge theory 

This appendix is intended to give a detailed description of AA = 2 quiver gauge theory, 
including details which the authors have not found elsewhere in the literature. 

A.l Relation to AA = 4 SYM theory 

In this section we give a detailed description of how N = 2 U{N)^ quiver gauge theory 
can be obtained by applying a Z,m projection to AA = 4 U{NM) SYM theory. 

Consider Type IIB string theory and introduce a stack of NM coincident D3-branes into 
the 10-dimensional (initially flat) spacetime. It is well known that the low-energy effective 
field theory of open strings with endpoints attached to the D3-branes is 4-dimensional 
AA = 4 SYM theory with gauge group U{NM). The space transverse to the world volume 
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of the D3-branes is = C'^ which has the isometry group SO (6). Now we consider the 
action of the subgroup Zm on given by 



The group Zm is called the orhifold group. We will denote the resulting quotient of by 
C^/Ztv/ where it is implied that the action of Zjv/ on is always that given in Eq. (A.l). 



Consider now open strings living on the stack of D3-branes where the transverse space 
is C'^/Zm- The low-energy effective field theory is no longer J\f = A U{NM) SYM theory. 
This is because associated with the orbifold group action (|A.l] ) on the coordinates of 
there is an orbifold group action on the scalar fields and their superpartners (to be defined 
below), and we must require that all quantum fields of = 4 SYM theory be invariant 
under this action. The gauge theory obtained from = 4 SYM theory by truncating the 
Hilbert spaces of quantum fields to Zjvf-invariant fields is called M = 2 quiver gauge theory. 

The i?-symmetry group of = 2 quiver gauge theory is SU{2)ji x U{l)ji. This is 
shown explicitly in Appendix A. 2 where the Lagrangian density of the quiver gauge theory 
is expressed in terms of SU{2)fi x U{1)r invariants. The quiver gauge theory thus indeed 
has M =2 supersymmetry. 

The orbifold group action ( |A.l] ) breaks the gauge group U (NM) of the = 4 theory 

into 

C/(iV)(i) X [/(iV)(2) X • • • X C/(iV) W (A.2) 

which is thus the gauge group of AA = 2 quiver gauge theory. We can see this as a manifes- 
tation of the fact that the quiver gauge theory is a low-energy effective field theory of open 
strings. Indeed, each of the M copies of C^/Zm embedded in will contain coincident 
D3-branes, and an open string can attach its endpoints to any of the stacks. Finally, to 
conclude the enumeration of the symmetries of M = 2 quiver gauge theory, we note that it 
is known to be a conformally invariant theory like the parent AT = 4 SYM theory p5| . 

In order to define the action of the orbifold group Zm on the AA = 4 SYM fields we 
first set up some notation. First Zm is embedded into U (NM) by defining the twist matrix 
7 = diag(l,u;, . . . ,u>'^^^^) and mapping Zm ^ k i — > G U{NM). (Note that the entries 
uj^ of 7 are really N x N matrices. )^^ The action of Zm on the AA = 4 SYM fields is then 

0^(7')^(P-'A)7' (A.3) 

where p ■ cp equals a phase times the field (p. For the scalar fields the phase is determined by 
their identifications with the zi,Z2 and z^ directions in and comparing with ( A.l| ). For 



the gauge field the phase is 1. For the spinor fields the phase equals that of their bosonic 
superpartner. Thus the condition for the AA = 4 SYM fields </) to be invariant under the 
action of Zm is 

^ = iHp-<P)1- (A.4) 



^^Note that this representation of Zm satisfies Tr 7*^ — for all k £'Lm \ {0}. As pointed out in Ref. [Q, 
this is needed for consistency (the cancellation of one-loop open string tadpole diagrams). 
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In the following we will obtain the Lagrangian density of = 2 U{N)^ quiver gauge 
theory by rewriting the J\f = A U {NM) SYM Lagrangian density and require that all the 
fields satisfy the Zjvf-invariance condition ( A.4 ). 

We now consider = 4 U{NM) SYM theory on M x S"^ where the radius of 5^ is 
denoted by R. The scalar fields will couple conformally to the curvature of the 5^ through a 
quadratic term in the action. In the decompactification limit R ^ oo this term will vanish. 
The action of = 4 U{NM) SYM theory on M x 5^ equipped with a metric of Euclidean 
signature reads 

S^=^ = I d^x Tr + \{D,ct>'){D,<i,') + ^i?" - -/ [</'\ [<t>\ </'^] 

+ ^VvTm^mV-p - + i/3p\75<^'')' V',]) (A.5) 

where F^y = d^Ay — d^A^ + ig[A^,Ay] and = ig[A^, ■ ]. The traces are taken over 
the gauge indices. The indices have the ranges /x, = 0, . . . , 3; i, j = 1, . . . , 6; p, g = 1, . . . , 4 
and k = 1, . . . ,3. Here 0* are six real scalar fields and ipp are four 4-component Majorana 
spinors. Moreover, 7^ are the 4-dimensional 4x4 gamma matrices and and (3^ are 4x4 
matrices satisfying the relations 

{a'=,a'} = -25^'l4, = -25'='l4, [a^/3']=0. (A.6) 

Explicit representations can be given as 

f ' "A . «^=r T 1 . T r 1 , (A.7) 






The bosonic part of the quiver action 

To put the action of = 4 SYM theory in a form suitable for performing the orbifold 
projection we now define three complex scalar fields 

The fields 0* are Hermitian (since they transform in the adjoint representation of the gauge 
group U{NM)), so by Hermitian conjugation of ( A. 9 ) we find 

A=i=(<Ai-^(A'), B = ^{^^-icl>^), l=-L(^5-i</,6). (A.IO) 

The scalar field part of the AA = 4 SYM Lagrangian density written in terms of these fields 
takes the form 

= Tr (^D^ADf^A + D^BD^B + Z)^$L)^¥ 

+ R-^ Tr (^AA + BB + + £^='^ + L^=^ (A. 11) 
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where the D and F terms are, respectively, 



'-D 



]^g^l^[[A,A] + [B,B] + [^,^ 



2g' Tr [A, B] [A, B] + [A, [A, + [B, ^ [B, ^ 



(A.12) 
(A.13) 



The scalar fields and B can be identified with the zi,Z2 and Z3 directions of the 
(because they are the Goldstone bosons associated with breaking the translational invari- 
ance in the directions transverse to the D3-branes), so we have the orbifold group action 
p : A, -B) 1-^ {(^,uj~^A,ujB)^ and the condition for these fields to be invariant under the 
ZM-transformation is then 



r^'^(^^ = (^^ j'^A'y = ujA, 'y'^Bj = LL! ^B . 



(A.14) 



One easily checks that these conditions are satisfied by splitting the NM x NM matrix 
fields of the M = 4 U (NM) SYM theory up into M x M block matrices whose entries are 
N X N matrices: 



/A 



A,. 



A 



H2 



V 



A 



( Ai,2 

^2,3 



B 



( 

^2,1 
-63,2 

V 



A 



\ 



A 



\ Am,i 



(A/-1),M 





(A.15) 



B 





M,{M-1) 







/ $1 



V 



(A.16) 



M I 



Here ^^j, (j+i), j and are N x N matrices (where i = 1, . . . , M and we identify 

i ~ i + M). Inserting the Zjv/-invariant forms of A^^, A, B and ^ given in Eqs. ( |A.15 )-( A.16 ) 
into ( [A.lip -( |A.13 ) the scalar field part of the = 2 quiver gauge theory Lagrangian density 
reads 

M ( 



i=l 



+ Tr 



+ TV 



-^scalar = X/ 1 ( + ~ 



X ( + igA^iB(^i+i)^i - 



+ ig[Af,i, $i] d^<^i + ig[Af,i, 



+ Tr A 
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+ Tr 



2g'^ Tr - 



2/Tr 



(A.17) 



Inserting the form of given in ( A.15|) into (|A.5|) , the gauge field part of the N = 2 
quiver gauge theory Lagrangian density reads 



1 ^ 

r — - \^ Tr 



(A.18) 



i=l 



where of course = dfj,A[, — d^A^^ + ig[A^^,Al]. 

The fermionic part of the quiver action 

The fermionic part of the = 4 SYM Lagrangian density reads 



'1 ^ 
,0 -1, 



where the gamma matrices are given by 



(A.19) 



7/. = 




75 


= 70717273 




(1, ia) , 




EE (1, -ia) 



(A.20) 
(A.21) 



and representations of a'^ and /?'^ are given in Eqs. ( |A.7|) and ( |A.8|) , respectively. The fields 
^pp, p = 1, ... ,4 are 4-component Majorana spinors which can be decomposed in terms of 
2-component Weyl spinors as follows 



(A.22) 



where a = 1, ... ,4 is the spinor index on ipp. The Majorana spinors are related to their 
conjugates through the Majorana condition 



where the Majorana conjugation matrix is C ■ 



(A.23) 



^ 
e"^, 



with ei2 = —621 = —1. 
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Combining Eqs. (|A.22D and (|A.20|) - (|A.21|) one finds 



(A.24) 



— — 

Here the operator is defined by XpD^,Xq = 2{XpDf,Xq - iDf,Xp)Xq}- 

It will be useful for exhibiting the -R-symmetry of the quiver gauge theory to express 
the fermionic Lagrangian density in terms of the following Weyl spinors 



XA = Ai 



XB 



A, 



= A3 



(A.25) 



Here xa,Xb,'4'i''P^ are the respective superpartners of A, B, A^,^. Note here that the bar 
used over the spinors in ( [A.25| ) is understood to mean the Hermitian conjugate whereas the 
bar over the Xp in ( A.24| ) denotes the usual conjugate of Weyl spinors. Explicitly, letting 
a = 1, 2 be the spinor index and letting m, n be the gauge indices, 



(A 



1 )a,mn 



iXATa. 



{XA)a,: 



and 



(Al)a,mn — {{^iTa^mn) — (^l)* 



(x 



A)a,mn 



(A.26) 



(A.27) 



and analogously for XB,ip^ and ip. In particular, note that all the Weyl spinors XA,XB,ip^ 
and ^p have undotted indices. 

Inserting the definitions ( [A.25 ) into the decomposition ( [A.24| ) we can write the kinetic 
part of the fermionic = 4 SYM Lagrangian density (|A.19| ) in the form 



C 



,Ar=4,kin 
ferm 



iTr [xA T^D^iXA + XB t^D^^XB + V' T^D^ip + ^p<s> T^Df^ipq, 



(A.28) 



In order to find the potential part of the fermionic M = 2 quiver gauge theory Lagrangian 
density we first rewrite the analogous part of the = 4 SYM Lagrangian density ( |A.19 ). 
By inserting the explicit forms of the a^,f3^ matrices given in Eqs. 



lJ) into ( |Al9D 

and then decomposing the 4-component Majorana spinors into 2-component Weyl spinors 
according to ( A.22 ) and finally making the substitutions ( A.25 ), the = 4 SYM theory 
result may be expressed as 

-Cferm = iT^r [XAr^,Df,XA+YBr^D^,XB +'ipT^Df,'ip + -ip.i>T^D^,'ip.i>j 

+ TV (xa([A V-*] - [B,lp]) +x^{[A,lp] + [B,i;^^ 

-^{[AX^] - [B,xa\) - V'<i.([^,Xa] + [B,x^]) 
+ Xa(P,V^] - [B,^]) + xb{[A,^] + [B,^]) 

- i>{[A, xb] - [B, xa]) - Xa] + [B, xb]) 

+ xa[^,xb] - xb[^,xa\ + V'[$,V'«i>] - V'^i^,?] 
+ xa['^,xb] - xb['^,xa] + 



(A.29) 
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The Weyl spinor fields XA, XSj V'"^) "0 are the respective superpartners of A, B,^, A^. There- 
fore they must satisfy the Zjvf-invariance conditions 



(A.30) 



One easily checks that these conditions are satisfied by splitting the NM x NM matrix 
fields of the Af = 4 U (NM) SYM theory up into M x M block matrices whose entries are 
N X N matrices: 

/ XA,i \ 



/V'l 



V'2 



V 



( 

XBA 



XB 



XA 



Xb,m 



XA,2 



XB,2 







V 



\XA,M 



XA,M-1 





(A.31) 



V 



(A.32) 



XB,M-l / 

Here '4'i,XA,i,XB,i and tp^^i are N x N matrices (where i = 1,...,M and we identify 
i ~ i + M). Inserting the Zj\,/-invariant forms of ip,XA,XB and ip^ given in Eqs. ( A.31 )- 
( |A3^ ) into (|AT29|) , the spinor field part of the Af = 2 quiver gauge theory Lagrangian 
density reads (summation over i = 1, . . . , M implied) 



^ferm = « Tr \^XA,i T^D^XA,i + XB,i T^D^,XB,i + ^l^i T^j,D^j,1pi + T^Df^Tp-S,,; 

+ -^Tr (^xI7^i,(i+i)V''i>,(i+i) - XXlV'*,i-4j^(i+i) - Xi;7-B(j+i),iV'i+i + Xi;iV'i^(i+i),: 



+ XB,i^j,(j+l)V'i - XB,i'lpi+lA,{i+l) + XB,i^(j+l),jV'<I>,j - XB,jV'3>,(i+l)^(j+l),i 

- + 1p<l>,{i+l)XA,iAi,{i+l) - + 

+ XA,j^i,{j+l)V''I',i - XA,iV'$,(i+l)^j,(i+l) - XA,i-B(i+i),iV« + XA,iV'j+l-B(j+i)^j 



-V'i^i,(i+l)XB,j + V'j+lXB,i^j,(i+l) + V'i+l-B(j+l),iXA,j - fpiXA,iB(^i+i)^i 



- '>P^,ii+l)A,{i+l)XA,i + ■0*,iXA,i^j,(j+l) - '0<I>,j^(i+l),jXB,i + '0*,(i+l)XB,i^(i+l),: 
+ XA;7^iXB7 - XA,iXB,i^j+l - XB^^i+lXAi + XB/iXA,i^i 
+ XA,i^i+lXB,i - XA,iXB,i^i - XB,i^iXA,i + XB,iXA,i^i+l 



(A.33) 



We conclude that the Lagrangian density of AA = 2 U{N) quiver gauge theory is 



C = C 



scalar 



;auge 



+ Ae 



(A.34) 



where jCscaiar ) -^gauge and >Cferm are given in Eqs. (|A.17| ), (|A.18| ) and ( |A.33| ), respectively. 
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A. 2 i?-symmetry 

The Lagrangian density of = 2 quiver gauge theory (given in Eqs. ( A. 3^ ), ( [A.17D , ( [A.18 ) 
and ( A. 33 )) is invariant under global 5C/(2)r x U{1)r transformations. The C/(1)_r factor 
of the i2-symmetry group acts on the fields as 

''^<^i (A.35) 

(A.36) 
(A.37) 

The U{1)r transformations of the Hermitian conjugate fields are obtained by flipping ( — > 
— C. The Lagrangian density is manifestly invariant under the U{l){i transformation. 

We now move to consider the SU{2)ji transformations. Define the 2-component spinors 









— ' , $i — > 




XA,i — 




XB,i — > e-'^/'^XB,i 




Ipi — 







(A.). ^ , (Xlr - t''^'^'' ■ (A.38) 

Under a G SU{2)r these spinors have the transformations 

{\i)a aj' {Xi)b (A.39) 
(A-)" iW^b^- (A.40) 

Note that (A^) 

a — ^abi^i)^ has the transformation 

(\)a eabaJ'e'''(X;)d = a^''{X)d (A.41) 

where the equality follows by using a G SU{2)ji. Thus, {Xi)a and {Xi)a are SU{2)ji doublets. 
To exhibit the SU{2)r invariance of the Lagrangian density we define SU{2)ii invariants 
such as 

(Ai)a (Ai)" = -e°'''{\i)a{Xi)b = - (A. 42) 

and write the Lagrangian density in terms of these. For Af = 2 quiver gauge theory the 
bifundamental scalars and the adjoint fermions are organized into SU{2)r doublets as 
follows 



The scalar field Lagrangian density written in terms of the SU (2) r doublets takes the 
following form^^ 



^'^Note that the term R ^ Tr [e''''{Xi)a{Xi)b + "l>i$i) describing the conformal coupling of the scalar fields 
to the curvature has been omitted here. 
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M 



1=1 



- 2g^ Tr [e''\Xi)aiX)b - e'^^(A~T)„(A,_i)6)' 

- 2g^ Tr (e"^(Ai)„(A^)fee'^<^(A~T)e(Ai_i)<i + e«^(A,_i)„(A,)6e^'^(Al)e(A~T)d) 
+ 2ff2 Tr [e''\XUXi)be''''(XlUXi)d + e''\XMX)be'''i\UX)d) 

- 2^2 TV (e''\Xi)a^i+i(Xi)bWi + e''\Xi)a^i(Xi)b^i) 

+ Tr (e"^(A0a^$i+i(A^)6 + e''\Xi)a(XlWi^z)] ■ (A.45) 

The spinor field Lagrangian density written in terms of the SU {2)ji doublets takes the 
following form 



i=l 

+^TV (e-^{XA,i(Al)e, (Xl)4 + e'=''{xA,i, ixr^UX)d} 



+ e''^{xi;7(A^)c, (x*+i)4 + ^^'^{xAi, {Xi)c{Xi)d} 
+ e"^{xB,^(A^)c, (XITT)4 + e'"^{xB,^, ixl)c{Xi)d} 

- ^'"^{XB^iiXDc, {Xi)d} - f'^^ixB^i, {Xi+l)c(Xi)d} 

+ e^{{x^)c^^, {x^)d} + e^'{fe)cll, (x7)4 
+ {XA,i^i+l, XB,i] + {xXI'S'n XB^} 

- {XB,i'^i, XA,i} - {xB^^i+1, XA~i}) 



(A.46) 



These results are conveniently summarized in Table A which lists the i?-charges of all the 
fields in A/" = 2 quiver gauge theory. 













XA,i 


XB,i 






U{1)r 








1 





1 

2 


1 

2 


1 

2 


1 

2 


SUi2)R 


1 

2 


1 

2 














1 
2 


1 
2 



Table A. R- charges for the bosonic and fermionic fields 

Here the generators of 5u(2)/j are taken in the fundamental representation and chosen as 
^(<7a;, <7y, (Tz). The i?-charges of the corresponding Hermitian conjugate fields are obtained 
by simply changing the signs of the U{l)ji and SU{2)r charges. 



-47- 



B. Bosonic and fermionic fluctuation determinants 



In this appendix we present some technical details of the computation of the 1-loop quantum 
effective action given in Section 5. More specifically, we explain here how to evaluate the 
fluctuation determinants arising from path integrating over the fluctuating fields. 



B.l Bosonic case 

The fluctuation operators D^", Dg" and D^" in Eqs. (|]2l|)-(|]2|) are given as below. 



and 



and 



^2 (^(a(i-l),i)rm (l(i-l),i)m)n + (^i,(i-l))nri (^i,(i-l))mm 



for j = i — 1 



+ (^i,('i-l))nn (&i,(i-l))rm + (^(i+l),i)mm 
+ {i4'i)nn - (0i)mm) (((?!'i)nn " (0i)mm) 



(B.l) 



-2 



((fli,(i+l))rm (ai,(i+l))mm + (^(i+l),i)nn (&(i+l),i)n"n) 



-2 (^(ai_(i-|_i))nn (a(i-l),i)mm + (fe(j+l),i)nn (6i,(i_l))i: 

+ 2 ^(ai^(i-|_i))rm (ai,(i+l))nn + (ai^(i+l))mm (ai,(i+l))mr 
+ (^(i+l),i)nn (6(i+l),i)jm + (fe(i+l),i)mm (fe(i+l),i)r 
+ ((0 i+1 )nn i4>i)mm) ((0j+l)nn " mr? 



for j = i 
for j = i + 1 

for j — i — 1 



(B.2) 



for 7 = i 



-2 (^(a(.i+l),(j+2))nn (ai,(i+l))mm + (&(i+2),(i+l) (^(i+l),i)mm) for j — i + 1 



-2 



((a(i-l),i)nn (ai,(i+l))mm + (&i,(i-l))rm (^(i+l),j)mm^ for j — I - 1 

92 _ 2z« - <+i)do + « - + 

+ 2 (ai,(j+l))riri + (^12,(1+1) )mm (fli,(i+l))mm 

+ - (</'i+l)mm) ((0i)n« " (0i+l)mm)) for j = i 

^2 (^(ai,(i+l))nn (a(i+l),(i+2))mm + (&(i+l),i)rm )mm) for j = i + 1 



(B.3) 
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and 



~2 ^(a(i-l),i)nri io-(i-l),i)mm + (^i,(i-l))rm 1) )mm^ 

+ 2 (^(ai,(i+l))nn (ai,(i+l))rm + (a(i-l),i)mm (a(i-l),i)r™ 



+ (^i,(i-l))nn (^i,(i-l))nn + (^(i+l),i)mm 
nn (0i)mm 



2 (^(ai,(i+l))nn )mm + (^(i+l),i)nn (^(i+l),i)mm^ 



for J = i — 1 



(B.4) 



for j ~ i 



for J = i + 1 



In the general vacuum ( |5.10| )-( |5?T3 ) these operators are tridiagonal, periodically continued 
matrices (assuming M > 3). The determinant of this class of matrices was considered in 
Ref. ||59| (Appendix B) who found the following result, valid for M > 3: 



detn' 



/ {r~\mn\ fr-\mn\ { r-\mn\ 

i=M \ ^ ^ 



+ (-l)*^+Hr n ^ q'^*"'^ 

i=M \ 







1\ 



(B.5) 



The inverse order of the initial and final indices on the product symbol indicates that the 
matrix with the highest index i is on the left of the product. 

Fortunately, in the vacuum ( |5.14| )-( |5?T^ ) the fluctuation determinants take a much 
simpler form. Namely, using (|5H)-(|5l^), the operators n™"^ n^"^ ng"^ nmri ^^j. g^^^ 
m, n) can be written in the particular form below, and there is a simple closed expression 
for the determinant.^^ That is, defining have the determinant formula 



det 



-UJ 



h 

UJ 7] 



_^fc(M-2)^ 



M 



(B.6) 



i=l 



Note in particular that the phases uj'^ on the left hand side cancel out. Therefore, for any 
value of A: E Z in (|]l|)-(|5l|), one obtains the same result for the fluctuation determinants. 



B.2 Fermionic case 

In order to compute the fluctuation determinant arising from path integrating over the 



^*To prove the formula, note first that the powers of o;*^ appearing in the super- and subdiagonal mutually 
cancel according to (B.5), so the determinant is independent of the value of k. Putting k — Q, the formula 



(B.C) is a special case of Eq. (A.l) in Ref. ||^ 
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fermionic fluctuations we must first introduce some notation: 



and 



/ 1 

oj^ 

■■■ 



1 

uj^ 



(B.7) 



wi = a = {A) , 
W2 = b = {B), 
^3 = , 



^^1 = 

^^3 = 1m 



(B.8) 
(B.9) 
(B.IO) 



where it is implied that A,B,^ take the Zm projection invariant forms given in ( A.15|) - 
(|Al^ )- 

The fluctuation operator D^" in (|]26|) is as given below (where c,d = 1, 2, 3) 



2dp,r,(d, + d^o(a, a, ))6^, ^^^^^^^^ _ _ ^^^^^^^^ 



D" 

■"-'7, 



(B.ll) 

The reason why the Wc entries labelled by the gauge index m have additional factors of 
compared to the entries labelled by n comes from the commutator structure of the Yukawa 
coupling (see ( A.29| )). Namely, when taking the trace over the gauge indices, the Wc entries 
labelled by n correspond to the terms where a scalar field appears between two spinor fields, 
whereas those labelled with m correspond to the terms where the scalar field appears to 
the right of both spinor fields. After substituting the orbifold projection invariant forms 
given in Eqs. ( |A.15D -( [Al6|) and (|A.31| )- (|A.32| ), the bifundamental scalar VEV's will couple 
different pairs of spinor fields depending on whether the VEV appears between the spinor 
fields or to the right of them in the Yukawa coupling. Since the scalar VEV's are mutually 
related through the vacuum ( ^.14D -( ^?T7| ), this can be compensated for by appropriately 
multiplying factors of 

To compute the result of the path integrations it is convenient to define (for a fixed c) 



Gc = ((^c)rin (^^c)rin) i^\^c)mmS^c ^\Wc)mmS^c ) • 



Noting that 



(B.12) 
(B.13) 

(B.14) 
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one finds, by using tlie (anti) commutation relations ( |A.6|) for and /J'^, that the result of 
the path integrations over the fermionic fluctuations (Ap)j, (Ap)j is 



\{\{a\a%rF,Fd + {a\i5%rF,G,, - . l5%rG cG a) ) (B.15) 



det(D^") = det ( - (9^ + i<5^o(«r " "D) 

= det [{id^ - («r - «r))' + i(^cFd - G,Gd) 5cdS, 
= det Aij . 

Here we have defined the M x M matrix (labelled by i, j = 1, . . . , M) 



(B.16) 
(B.17) 



-2 ^(ai_2)rm (ai,2)mm + (fe2,l)nn (&24)mm^ ^'''^ 

+ 2 ( (ai^2)nn (ai,2)nn + {cLia)mm (^1,2) mm 

+ (^2,l)nn (^2,l)nn + (^2,l)mm (^2,l) mm 

+ ((0l)nn - ((/>l)mm) {i4>l)nn " (0l) mm 7 I 



for j = i — 1 



(B.18) 



— 2 ^(ai^2)nn (ai,2)mm + (i'2a)nn (fe2a)mm^ (^^^ 



for J = i 
for j — i + 1 



where we have used ( |5.14| )-( |5?T^ ) to arrive at the equality ( B.17| ). Applying the determinant 
formula ([B.6| ) and using ( 5.1^ )- (|5.17 ) again one finds, after taking the traces over the 
fermionic Matsubara frequencies = (^^+^)^ and over the spherical harmonics, the 
expression (|5.27 ). 
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